Real-time renormalization group in frequency space: A two-loop analysis of the 
nonequilibrium anisotropic Kondo model at finite magnetic field 



Herbert Schoeller^'^ and Frank Reininghaus^' ^ 

' Institut fiir Theoretische Physik, Lehrstuhl A, RWTH Aachen, 52056 Aachen, Germany 
JAR A- Fundamentals of Future Information Technology 
(Dated: July 23, 2009) 

We apply a recently developed nonequilibrium real-time renormalization group (RG) method 
in frequency space to describe nonlinear quantum transport through a small fermionic quantum 
system coupled weakly to several reservoirs via spin and/or orbital fluctuations. Within a weak- 
coupling two-loop analysis, we derive analytic formulas for the nonlinear conductance and the kernel 
determining the time evolution of the reduced density matrix. A consistent formalism is presented 
how the RG flow is cut off by relaxation and dephasing rates. We apply the general formalism 
to the nonequilibrium anisotropic Kondo model at finite magnetic field. We consider the weak- 
coupling regime, where the maximum of voltage and bare magnetic field is larger than the Kondo 
temperature. In this regime, we calculate the nonlinear conductance, the magnetic susceptibility, 
the renormalized spin relaxation and dephasing rates, and the renormalized g factor. All quantities 
are considered up to the first logarithmic correction beyond leading order at resonance. Up to a 
redefinition of the Kondo temperature, we confirm previous results for the conductance and the 
magnetic susceptibility in the isotropic case. In addition, we present a consistent calculation of the 
resonant line shapes, including the determination whether the spin relaxation or dephasing rate cuts 
off the logarithmic divergence. Furthermore, we calculate quantities characterizing the exponential 
decay of the time evolution of the magnetization. In contrast to the conductance, we find that the 
derivative of the spin relaxation (dephasing) rate with respect to the magnetic field is logarithmically 
enhanced (suppressed) for voltages smaller (larger) than the renormalized magnetic field, and that 
the logarithmic divergence is cut off by the opposite rate. The renormalized g factor is predicted to 
show a symmetric logarithmic suppression at resonance, which is cut off by the spin relaxation rate. 
We propose a three-terminal setup to measure the suppression at resonance. For all quantities, we 
analyze also the anisotropic case and find additional nonequilibrium effects at resonance. 

PACS numbers: 05.10.Cc, 72.10.Bg, 73.63.Nm 



I. INTRODUCTION 

One of the basic unsolved problems of dissipative quan- 
tum mechanics and quantum transport through meso- 
scopic systems is the nonequilibrium Kondo model. In 
its simplest version, it consists of a spin-i system cou- 
pled via exchange processes to the spins of two fermionic 
reservoirs, which are kept at two different chemical po- 
tentials ^L/fl — ±6^, see Fig.[T] Besides the importance 
of the Kondo model for many aspects of strongly cor- 
related Fermion systems (see Ref. |l| for an overview), it 
was suggested to realize this model in transport experi- 
ments through quantum dots.^ This has been achieved^ 
with the particular advantage of full control over all pa- 
rameters such as temperature, voltage, magnetic field, 
and exchange couplings. The central idea is to lower 
temperature and bias voltage such that only one single- 
particle level of a quantum dot will contribute to trans- 
port. Adjusting the gate voltage such that charge fluc- 
tuations of this level are suppressed (Coulomb blockade 
regime), the dot can either be occupied by a spin up 
or a spin down electron, and the spin can fluctuate via 
second-order cotunneling processes, leading precisely to 
the model depicted in Fig. [TJ In this realization, one ob- 
tains an antiferromagnetic exchange coupling of the order 
J '-^ pt^ /Ec (in dimensionless units), where p is the den- 
sity of states in the reservoirs, t denotes the tunneling 
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FIG. 1: A spin-i quantum system coupled via exchange to 
two reservoirs. Jl = Jll and Jr = Jrr involve exchange 
between the spins of the left/right reservoir with the local 
spin, and Jlr = Jrl transfers a particle from one reservoir 
to the other during the exchange process. 



amplitude (hopping parameter) between the leads and 
the dot, and Ec is the energy of the virtual intermedi- 
ate state (charging energy). In equilibrium, the Kondo 
model has been analyzed by various many-body methods 
(for an overview, see Ref. [l|) and can be solved exactly 
by Bethe ansatz"* or conformal field theory.^ Powerful 
numerical techniques such as the numerical renormaliza- 
tion group have been developed^ii from which all ther- 
modynamic and spectral properties can be calculated. 
The basic physics can already be understood from poor 
man scaling arguments^ If all exchange couplings are the 
same Jl = Jr = Jlr = Jo, one obtains the renormaliza- 
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tion group (RG) equation 
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Here, Tk is the Kondo temperature and J denotes an 
effective exchange coupling corresponding to an effective 
band width A of the reservoirs {D ~ Ec is the origi- 
nal band width). For the antiferromagnetic case J > 0, 
one obtains an enhancement of the exchange coupling by 
reducing the effective band width until A reaches Tk, 
where J diverges. This indicates a logarithmic enhance- 
ment of the linear conductance for temperatures above 
the Kondo temperature (weak-coupling regime), and a 
complete screening of the dot spin by the reservoir spins 
below Tk (strong-coupling regime). In the latter case, 
it can be shown that the conductance becomes unitary, 
i.e., G = 2e^//i, at zero temperature T = and zero bias 
voltage V = Q. 

In nonequilibrium, the Kondo model is not yet solved 
completely. Some numerical techniques have already 
been developed to describe either the time evolution of an 
out-of-equilibrium initial state by using time-dependent 
numerical renormalization group (TD-NRG) (Ref.i) or 
to describe the stationary current in the presence of a 
bias by using NRG in a scattering wave basisi^ or using 
an iterative real-time path integral approach. ^-"^ Scatter- 
ing wave Bethe ansatz methods based on the Lippmann- 
Schwinger equation are currently under way to solve 
the single-impurity Anderson model in noncquilibriumjii^ 
Concerning analytical RG methods, it has been empha- 
sized that it is important to understand how the RG 
equations are cut off by the bias voltage and by relax- 
ation and dephasing rateS i^'^'^^i^^ The problem can be 
stated as follows. Performing perturbation theory in the 
bare exchange coupling Jq (disregarding for the moment 
the differences between J^, J^j, and Jlr), logarithmic 
terms occur, which, at zero temperature and depending 
on the physical observable under consideration, have the 
form [we indicate only the order in the exchange coupling 
leaving out numbers of 0(1) and other physical energy 
scales from V and /iq in the prefactor] 
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with k>l>0, n = 0, ±1, ±2, . . . and to = 0, ±1. Here, 
/ is the current, T1/2 are the spin relaxation/dephasing 
rates, g is the renormalized g factor, and M denotes 
the magnetization. D is the band width of the reser- 
voirs and /iQ is the bare magnetic field, n denotes 
the number of particles transferred between the reser- 
voirs and TO = 0,±1 characterizes whether spin flip pro- 
cesses occur or not. The points nV = mho correspond 



to resonance positions, where certain higher-order pro- 
cesses are allowed by energy conservation, e.g., the value 
n = TO = 1 corresponds to the onset of inelastic cotun- 
neling processesJ^ Higher-order terms with n > 2 have 
so far not been discussed in the literature for the Kondo 
model, but are generically expected'^^ (for other mod- 
els in the charge fluctuation regime, corresponding terms 
have been calculated in Ref. 18). For given perturbation 
order k, the allowed values for / depend on the value of 
n. For n = 0, 1 it is known that I = 0,1, ... ,k. Even in 
the weak-coupling regime 
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the logarithmic terms can lead to a breakdown of per- 
turbation theory if Jo In ^ is small enough such that 

JqIu;!^ — ^ — JqIu^ ^ 0(1). Furthermore, at reso- 
nance nV — mho, the logarithmic terms even diverge. 
Therefore, it is necessary to resum the logarithmic terms 
in an appropriate way using RG and, at the same time, 
introduce the physics of relaxation and dephasing rates 
to cut off the divergencies at resonance. This idea has 
been proposed in Ref. 19. In a first step, within a stan- 
dard poor man scaling approach, one resums all leading 
order logarithmic terms of the form 
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where r = 0,1,2 depends on the physical observable un- 
der consideration [see Eq. ([3])]. This leads to an effec- 
tive exchange coupling Jc, given by Eq. ^ evaluated at 
A = Ac: 
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In a second step, one tries to expand the physical observ- 
able systematically in powers of the effective coupling 
constant Jc, leading to a series with terms similiar to 
Eq. ([3]), but with the replacements 13 — > Ac and Jo — + Jc. 
If, in addition, one cuts off the resonances by an appro- 
priate relaxation or dephasing rate F, a new series of the 
form 
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is expected, where h is the renormalized magnetic field. 
As pointed out in Ref. [l^ this perturbation series in Jc 
is well-defined for Jc ^ 1, because the maximum value 
of the logarithm at resonance is given by In ^ In Jc, 
where we have used the rough estimate 

f^jfAc (8) 

based on a simple dimensional analysis. Therefore, at 
resonance, we expect terms of the form 
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which are not dangerous since In Jc <C 1 if Jc <C 1 
(note that k > I). The leading order result is the term ^ 
(denoted as one- loop in this paper), whereas the first 
logarithmic correction ~ Jc~^^ In Jc is the first subleading 
term (denoted as two- loop in this paper). 

The purpose of the present paper is to present a well- 
defined two-loop RG approach to work out the above- 
described procedure. Thereby we will apply a recently 
proposed real-time renormalization group method in fre- 
quency space (RTRG-FS).^^ This method has the ad- 
vantage that formally exact RG equations can be set 
up in nonequilibrium which include the relaxation and 
dephasing rates in all resolvents appearing on the right 
hand side (r.h.s.) of the RG equations. Furthermore, 
in leading order, precisely poor man scaling equation ([!]) 
is obtained. This provides the possibility to proceed in 
two steps: First one expands the exact RG equations 
systematically around the poor man scaling solution for 
D > A > Ac, and, in the second step, one solves the 
RG equations perturbatively in Jc for Ac > A > 0. As 
we will show in this paper, in both steps two-loop terms 
are important to obtain the first logarithmic corrections 
beyond leading order. In the first step, two-loop terms 
arising from higher-order terms on the r.h.s. of the RG 
equation generate terms ~ Jc^^ 7^ which depend only 
weakly on the voltage and are incorporated into a re- 
definition of the coupling constant Jc (or, equivalently, 
the Kondo temperature). These terms lead to an over- 
all increase (or decrease) in the physical observable un- 
der consideration but show no interesting dependence on 
voltage or magnetic field. In contrast, in the second step, 
logarithmic contributions of form ([7]) are generated which 
give rise to a logarithmic enhancement (suppression) at 
resonance. This means that the various two-loop terms 
(leading all to terms of the same order of magnitude) can 
be systematically divided into important and unimpor- 
tant terms concerning their dependence on voltage and 
magnetic field. In accordance with results of Refs. [H 
and [20l . we emphasize that it is very important to in- 
clude the frequency dependence of the vertices generated 
from the first step D > A > Ac since this infiuences the 
prefactor of the logarithmic contributions calculated in 
the second step Ac > A > 0. 

In one-loop order (but including certain two-loop terms 
from the frequency dependence of the vertices), the con- 
ductance and the magnetization have been calculated 
previously for the nonequilibrium Kondo model. Pio- 
neering works are Refs. [la andll^, where the slave parti- 
cle approach was used in combination with the Keldysh 
formalism and quantum Boltzmann equations. The RG 
was formulated only on one part of the Keldysh contour 
and a real-frequency cutoff was used. In these works, it 
was investigated how the voltage and the magnetic field 
cut off the RG flow and it was emphasized that the fre- 
quency dependence of the vertices is important to obtain 
the first logarithmic contributions beyond leading order. 
In fact, the result of the present paper concerning the 
conductance and the magnetization is precisely the same 



as that of Refs. [l^ and [13, up to the redefinition of the 
Kondo temperature. This means that we will prove in 
this pap er that all two-loop contributions neglected in 
Refs. 15 and '20 do not influence the prefactor of the first 
logarithmic correction beyond leading order. Further- 
more, in Refs. [isl and [20l. diagrams connecting the upper 
with the lower part of the Keldysh contour have been ne- 
glected. Within these works, it was therefore not possible 
to describe the cutoff of the logarithmic terms by relax- 
ation and dephasing rates on a full microscopic level. In 
this paper, we will show how this can be achieved within 
RTRG-FS, which provides a consistent formalism to cal- 
culate the line shape at resonance and to see whether 
the spin relaxation rate Fi or the spin dephasing^ rate 
r2 cuts off the logarithmic divergencies (in Ref. 29, the 
latter question was adressed only in bare perturbation 
theory and zero magnetic fleld). In addition, in this pa- 
per we will also discuss the anisotropic case and calcu- 
late the spin relaxation/dephasing rates together with 
the renormalized g factor up to the flrst logarithmic con- 
tribution (corresponding to a three-loop calculation in 
conventional RG methods). Besides the known reduc- 
tion in the magnetic field in first order in J^ we find 
that the renormalized magnetic field in second order in 
J is proportional to logarithmic terms similiar to Eq. ([7]) 
with a significant dependence on voltage and magnetic 
field. We propose an experimental setup with a weakly 
coupled third lead to measure the voltage dependence of 
the renormalized g factor. Moreover, we find that the log- 
arithmic terms of r2 and g are controlled by Fi, whereas 
those of Fi are controlled by F2. In the anisotropic case 
this leads to the effect that the logarithmic resonances 
of become sharper with decreasing J-*- since Fi does 
not contain any terms proportional to J^ in second or- 
der. Furthermore, we will show that the susceptibility de- 
pends only weakly on the tranverse coupling and the 
logarithmic resonances even survive in the limit J^ — > 0. 
The anisotropic Kondo model has recently been proposed 
to be realizable in low-temperature transport through 
single molecular magnets^ ^ and experiments are starting 
to investigate such systems. Since the transverse cou- 
pling is induced by small magnetic quantum tunneling 
terms, giving rise to rather small Kondo temperatures, 
the susceptibility might be an interesting physical quan- 
tity to measure signatures of the Kondo effect even for 
very small values of J"*". 

Using flow equation methods, ^'^ a consistent two- 
loop approach including the cutoff by spin relax- 
ation/dephasing rates has been presented in Ref. [2^ for 
the isotropic Kondo model in the absence of a magnetic 
field. Within this method, the cutoff from the rate F 
occurs due to a competition of certain one-loop and two- 
loop terms on the r.h.s. of the RG equation for the ver- 
tex. This is a completely different picture compared to 
RTRG-FS, where the cutoff parameter F together with 
the voltage occurs already in the one-loop terms as an ad- 
ditional term in the denominator of the resolvents. Thus, 
the RTRG-FS method is closer to conventional poor man 
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scaling and the physics of relaxation and dephasing rates 
occurs naturally as a resummation of a geometric series 
similiar to self-energy insertions in Green's function tech- 
niques. In this sense the RTRG-FS method proves that 
conventional scaling equations (properly generalized to 
the Keldysh contour) can account for the cutoff of the 
RG flow by rates and the voltage. Furthermore, the 
RTRG-FS method provides a generic proof in all orders 
of perturbation theory in the renormalized vertices that 
the various cutoff scales are the physical relaxation 
and dephasing rates governing the time evolution of the 
reduced density matrix of the quantum system. 

The RTRG-FS method used in this paper has been 
proposed in Ref. [l3 and is a natural generalization of 
an earlier developed real-time RG method.— i^^i^l As de- 
scribed in detail in Ref. [Tt! . many technical improvements 
have been incorporated, the main ones being a formula- 
tion of the RG in pure frequency space, integrating out 
the symmetric part of the Fermi distribution function be- 
fore starting the RG, and formulating the nonequilibrium 
RG on the imaginary frequency axis. As a consequence, 
the rates determining the cutoff of the RG flow obtain 
the right scale, and it is possible to show generically that 
relaxation and dephasing rates cut off the RG flow in all 
orders of perturbation theory and within all truncation 
schemes. Furthermore, the dependence on the Keldysh 
indices can be completely avoided, and one can calcu- 
late the time evolution and the nonequilibrium station- 
ary state in pure Matsubara space without the need of 
any analytic continuation. The latter idea has first been 
proposed in Ref. [1^ in the context of nonequilibrium func- 
tional renormalization group within the Keldysh formal- 
ism. A particular advantage of the RTRG-FS method is 
that relaxation and dephasing rates occur naturally as 
the negative imaginary part of the eigenvalue of the ker- 
nel determining the kinetic equation of the reduced den- 
sity matrix, and do not arise from more involved combi- 
nations of self-energies and vertex corrections as in slave 
particle formalism.^® Furthermore, it is straightforward 
to calculate the time evolution from RTRG-FS since the 
RG gives directly the result for the kernel in Laplace 
space. 

For generic problems with spin and/or orbital fluctu- 
ations, it was described in Ref. [13 how to solve the RG 
equations analytically in the weak coupling regime up to 
one-loop order. In this paper we will provide the techni- 
cally much more involved two-loop case, which is neces- 
sary to calculate consistently the important logarithmic 
terms at resonance discussed above, see Eq. Q. The 
result will be applied to the calculation of the conduc- 
tance and the magnetic susceptibility of the anisotropic 
Kondo model in the presence of a magnetic field. Fur- 
thermore, we will also analyze quantities characterizing 
the time evolution of the Kondo model. Thereby, we will 
concentrate on the calculation of the dominant exponen- 
tial decay of the magnetization in two-loop order, which 
is determined by the spin relaxation and dephasing rates 
ri^2 and the renormalized magnetic field h. Finally, in 



addition to Ref. [13, we will generically show that pre- 
cisely these physical quantities control all resonant line 
shapes (in Ref. [13, the question whether the cutoff scales 
of the logarithmic terms are exactly identical to the phys- 
ical relaxation/dephasing rates was still open). 

The paper is organized as follows. In Sec. Ill A[ we will 
set up the generic model and the perturbative series. Sec- 
tion [TTB] is devoted to the general RG formalism and the 
derivation of the two-loop RG equations for an arbitrary 
quantum dot coupled via spin and/or orbital fluctuations 
to reservoirs. The systematic way to analytically solve 
these RG equations up to two-loop order is presented 
in Sec. Ill CI The general formalism is applied to the 
nonequilibrium Kondo model in Sec. lIIII In Sec. lIII A! we 
will set up the algebra in Liouville space needed to eval- 
uate all expressions explicitly, and Sec. IIIIBI describes 
the evaluation of the general two-loop equations for the 
Kondo model. The final results for the conductance, the 
magnetic susceptibility, the spin relaxation and dephas- 
ing rate, and the renormalized g factor are presented for 
the isotropic case in Sees. IIV M IIV Bl and IIV CI whereas 
the anisotropic case is discussed in section lTVDI We close 
with a summary in Sec. [Vj A list of all symbols used in 
this papers is presented in Sec. IVII 

We emphasize that a reader, who is not interested in 
the formal derivation but only in the final physical re- 
sults of the paper, can skip the following formal Sees. [Tlj 
and mil and can directly move over to Sec. IIVI where 
the notations needed to understand the results are again 
repeated. 

II. GENERIC CASE 

In this section, we describe a generic quantum dot cou- 
pled via spin and/or orbital fluctuations to several reser- 
voirs. In Sees. Ill Al and HTbI we introduce the basic nota- 
tions, the perturbative series and summarize shortly the 
setupof the RG equations as explained in more detail in 
Ref. [13. In Sec. Ill C[ we present a systematic way how to 
solve the RG equations analytically up to two-loop order 
in the weak-cougling regime (the one-loop case has been 
treated in Ref. [13) • Throughout this paper, we use units 
e = h = ks = ^■ 

A. Model and perturbative series 

Model. We consider a quantum dot with fixed charge 
in the Coulomb blockade regime where only spin and/or 
orbital fluctuations are possible via the coupling to exter- 
nal reservoirs. As shown in detail in Ref. ,27i , a standard 
Schrieffer- Wolff transformation leads to a Hamiltonian of 
the form 

H = Hres + Hs + V = Ho + V , (10) 

where Hres is the reservoir part, Hs characterizes the iso- 
lated quantum dot, and V describes the coupling between 
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reservoirs and quantum dot. They are given explicitly by 



Hs = 
V = 



j dw {uj + ^ia) a+^(uj)a-^{uj) , (11) 

E^«l^)(^l ' (12) 

S 
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Here, the creation (?? = +) and annihilation 

(77 = —) operators of the reservoirs, and u is an in- 
dex characterizing all quantum numbers of the reservoir 
states. In the absence of further symmetries, u contains 
the reservoir index a and the spin quantum number a 
(for two reservoirs and spin-i particles, we use the nota- 
tion a = L, R = zL and a =t, i= ±). w is the energy of 
the reservoir state relative to the chemical potential fia 
of reservoir a. The eigenstates and eigenenergies of the 
isolated quantum dot are denoted by \s) and Eg. The in- 
teraction V is quadratic in the reservoir field operators, 
which arises from second-order processes of one electron 
hopping off and on the quantum dot coherently (for neg- 
ative charging energies, also two electrons can hop off or 
on the dot^°). This keeps the charge fixed and allows 
only spin and orbital fluctuations. The coupling vertex 
9iiu.ri'u'{^i^') is an arbitrary operator acting on the dot 
states. It is written in its most general form, depend- 
ing on the quantum numbers and energies of the reser- 
voir states in an arbitrary way. However, as explained in 
Ref. [13, the RG approach can be set up in its most conve- 
nient form if one assumes that the frequency dependence 
of the initial vertices is rather weak and varies on the 
scale of the band width D of the reservoirs. Therefore, 
we will assume this in the following and introduce below 
[see Eq. a convenient cutoff function into the free 

reservoir Green's functions. 

To achieve a more compact notation for all indices, we 
write 1 = "qvLo and sum (integrate) implicitly over all 
indices and frequencies. The interaction is then written 
in the compact form 
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: • • • : denotes normal ordering of the reservoir field op- 
erators, meaning that no contraction is allowed between 
reservoir field operators within the normal ordering. A 
contraction is defined with respect to a grand-canonical 
distribution of the reservoirs, given by 
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Uiuj) = (e"/^° -I- l)-i = 1 - Ui-Lo) is the Fermi 
distribution function corresponding to temperature Ta 
(note that the chemical potential does not enter this for- 
mula since u; is measured relative to /Xq,). Furthermore, 
S{uj - Lu') is the S function in compact no- 
tation, and 1 = —f], V, ijj. The cutoff by the band width D 



can be introduced in many different ways into the reser- 
voir contraction. We use a Lorentzian cutoff and replace 
the contraction by 



ai av 611, p{uj) fairju;) 



with 



D2 



(16) 



(17) 



Within the normal ordering of Eq. (|14p . the field op- 
erators can be arranged in an arbitrary way (up to a 
fermionic sign) , therefore the coupling vertex can always 
be chosen such that antisymmetry holds. 
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Furthermore, due to the hermiticity of V, the vertex has 
the property 
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The particle current operator flowing from reservoir 
7 to the quantum dot is defined by H — -^-^N^es — 
—i[H, N^^g], where A^^es is the particle number in reser- 
voir 7. Using Eqs. pO|) and (|13p . a straightforward cal- 
culation leads to 



with 
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We are interested in the time evolution of the reduced 
density matrix ps{t) of the quantum dot and in the aver- 
age (/''') (i) of the current operator. Formally, they follow 
from the solution of the von Neumann equation 

:e-''^^'''°'> Ps{to)Pres , 
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{P){t) = TrsTvresi-iLi.)e-'''^'-*°^ psito)pre. m 
where 



L = [H,-l 



(25) 



are operators in Liouville space acting on usual opera- 
tors in Hilbert space via the (anti)commutator [A, B\± = 
AB ± BA. Initially, we have assumed that the density 
matrix is a product of an arbitrary dot part ps(to) and 
a grandcanonical distribution prea for the reservoirs. It 
is convenient to introduce the Laplace transform 
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FIG. 2: Example of a diagram for the reduced density matrix 
of the dot. The time direction is to the left. The dots rep- 
resent the interaction vertex G between the local quantum 
system and the reservoirs. Each vertex corresponds to two 
adjacent dots indicating that two reservoir field operators are 
associated with each vertex. The horizontal lines connecting 
the vertices denote the free time propagation of the quantum 



system, leading to the resolvents 



E+Xi 



in Laplace space. 



The lines connecting the dots are the reservoir contractions 
arising from the application of Wick's theorem. The dashed 
vertical lines between the vertices are auxiliary lines to deter- 
mine the energy argument Xi of the resolvents. 



Perturbative expansion. The next step is to expand 
expressions (|26)) and ((27| in the interacting part Ly = 
[V, •]_ of the LiouviUian and to integrate out the reservoir 
part. As outhned in detail in Ref. [ij, this leads to a dia- 
grammatic representation in Liouville space. We shortly 
summarize this procedure here. First, Ly can be written 
in the form 



(28) 



where Jf is a quantum field superoperator in Liouville 
space for the reservoirs, defined by {A is an arbitrary 
reservoir operator) 
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p = ± is the Keldysh index indicating whether the field 
operator is acting on the upper or the lower part of the 
Keldysh contour. G^, is a superoperator acting in Liou- 
ville space of the quantum dot, and is defined by {A is 
an arbitrary operator of the quantum dot) 
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Inserting the form (^5]) into Eqs. 1^5)) and (P7)) . expanding 
in Ly, and shifting all reservoir field superoperators Jf to 
the right, one can show that each term of perturbation 
theory can be written as a product of a dot part and 
an average over a sequence of field superoperators of the 
reservoirs with respect to Pres • 

Evaluating the latter with 
the help of Wick's theorem, one can represent each term 
of the Wick decomposition by a diagram, see, e.g.. Fig. [2] 
describing a certain process for the time evolution of the 
reduced density matrix of the dot. Each process consists 
of a sequence of interaction vertices G^f , between the dot 
and the reservoirs, and a free time propagation of the 
dot in between (leading to resolvents in Laplace space). 
Since the reservoirs have been integrated out, the vertices 
are connected by reservoir contractions (the solid lines 



without arrows in Fig. [2]). This means that the various 
diagrams represent terms for the effective time evolution 
of the dot in the presence of dissipative reservoirs. Each 
diagram for the reduced density matrix has the form 
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E- L 



(0) 



G 



E + X,-Lf 



■G ...G- 



where 



E + Xr- 



(0) 



-(0) 



[Hs 



G- 



E-L 



(0) 



Ps(io) , (31) 



(32) 



G = G^j^' indicates an interaction vertex, and 7 = jfj^^ 



is a contraction between the reservoir field superopera- 
tors, defined by 



hi' 



JP JP — J Tr IP JP n 

Sir p{uj)p' faivP^) ■ 



(33) 



To factorize the Wick decomposition, a fcrmionic sign 
has to be assigned to each permutation of reservoir field 
superoperators, indicated by the sign factor (—1)^" in 



Eq. ([3T|) . For each pair of vertices connected by two 
reservoir lines, a combinatorical factor ^ occurs, leading 
to the prefactor ^ in Eq. (|3T|) . The value of the frequen- 
cies Xi in the resolvents between the interaction vertices 
is determined by the sum over all variables x — r]{Lu + pa) 
of those indices belonging to the reservoir lines which are 
crossed by a vertical line at the position of the resolvent 
(see the dashed lines in Fig. [2|) . Thereby, the index of the 
left vertex has to be taken of the corresponding reservoir 
line, e.g., the diagram of Fig. [2] is given by (the obvious 
dependence on the Keldysh indices has been omitted for 
simphcity, i.e., 7^^- = 7; 



PiPj 



and G., ^G^^O 



E-L 



(0) 



716 723745 G12 U12 G34 G56 



-(0) 



2 77,10789 IIts Gg^io 



E-L 



(0) 



Psito) , (34) 



where the resolvents are defined by 

1 



ni...„ = 



with 



El. 



El. 



(0) 



Pz 



E 

i=l 



(35) 

(36) 
(37) 



As can be seen from example p4p . each diagram consists 
of a sequence of irreducible blocks (where a vertical line 
always cuts at least one reservoir line) and free resolvents 
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1/{E — L^g^) in between. Similiar to Dyson equations one 
can formally resum this series with the result 



Ps{E) 



E 



-eff 



Ps{h) 



with 



(38) 



(39) 



where the kernel S(-E) contains the sum over all irre- 
ducible diagrams. A similiar procedure can be used to 
calculate the average (P7|) of the current operator with 
the result 



{D){E) = Trs^^iE) 



1 



E 



reff 



(E) 



-iTTs^'r{E)psiE) , 



(40) 



where the current kernel E^(£') is defined similiarly to 
I](i?), but the first vertex G is replaced by the current 
vertex defined by 



(Tl\PP — A , nrPP 

such that, in analogy to Eq. 



(41) 



(42) 



Using Eq. ([3T|) . a certain diagram for the kernels has 
to be translated according to 



G 

r 



E + X,^Lf 



G ...G 



n 



1 



G , (43) 



where the subindex irr indicates that only irreducible 
diagrams are allowed where any vertical line between the 
vertices cuts through at least one reservoir contraction. 

The stationary solutions for the reduced density ma- 
trix and the current follow from the Laplace trans- 
form by = \imE^Q+{~iE)ps{E) and = 
\m\E->o+{—iE){D){E)^ and can be calculated from 



-/^(»0+)p|* 



(44) 



-zTrsS^(0+)pJ* . (45) 



In addition to previous formulations^'' of the pertur- 
bation series, we note that the diagrammatic series can 
be partially resummed by taking all closed subdiagrams 
between two fixed vertices together which contain only 
contractions connecting vertices between the two fixed 
ones. This has the effect that the resolvents in Eq. (l43l) 
are replaced by 



1 



E + X, 



(0) 



1 



E + Xi- Ly^{E + Xi) 



(46) 



i.e., the full effective Liouville operator occurs in the de- 
nominator. This means that Eqs. ([39|) and (|43|) turn into 
self-consistent equations for Lg^-^ [E) for any approxima- 
tion. Of course, the number of diagrams is reduced in this 
formulation. No diagrams are allowed anymore which 
contain closed subdiagrams between two vertices. 

When calculating diagrams with the replacement (|46|1 . 
one faces the problem that the frequency integrations 
cannot be performed analytically because the energy 
dependence of the effective Liouvillian is not known. 
This would require the solution of a complicated self- 
consistent integral equation. To avoid this, it is useful 
to formulate an appropriate approximation for the re- 
solvents which can be improved systematically. To de- 
fine this approximation, we write the resolvents in terms 
of the eigenvectors and eigenvalues of the Liouvillian 



n(z) 



where the projectors are defined by 

Pi{z) = \xi{z)){xi{z)\ , 



P^iz) (47) 



(48) 



and \xi{z)) and {xi{z) \ are the right and left eigenvectors 
ofr/^(z). 



L%ff{z)\x,{z)) = h{z)\x,{z)) , 



{xdz)\Llff(z) = K{z){x,{z)\ , 



(49) 
(50) 



with eigenvalues Ai(z). Assuming that \xi{z)) and {xi{z)\ 
have no poles (or poles with very large negative imagi- 
nary part so that they influence only the short-time be- 
haviour), the poles zt of the resolvent follow from the 
self-consistent equation 



(51) 



for all values of i. Expanding \i{z), \xi{z)), and {xi{z)\ 
around z = Zi, we see that the nonanalytic part of the 
resolvent is given by 



n(z) 



P^(^^) 



with residua (also called Z factors) given by 

1 



(52) 



(53) 



Equation (j52p defines our approximation which is the ap- 
propriate one to describe especially line shapes at reso- 
nance (analytic parts are expected to have no special 
features at resonance and will only lead to an overall 
perturbative shift of the background) . To avoid the sum- 
mation index i, we will write the approximation in the 
more compact form 



n(z) 



Ls 



(54) 
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where we use the convention that any function of Z and 
L5 is interpreted as 



(55) 



The eigenvalues Zi can be decomposed into real and imag- 
inary parts 



= hi 



r,; > 



(56) 



and are the poles of the original full resolvent n(z). Ac- 
cording to Eq. ([38l) . this resolvent describes the reduced 
density matrix in Laplace space. Therefore, the resolvent 
must be analytic in the upper half of the complex plane 
since otherwise solutions would exist in time space which 
are exponentially increasing and no stationary state can 
be reached. Thus, the negative imaginary parts Ti must 
be strictly positive and describe the various relaxation 
and dephasing rates of the different modes described by 
the eigenvectors \xi[zi)). Correspondingly, the real parts 
hi describe the oscillation frequencies of the modes, e.g., 
the effective magnetic field for the Kondo problem. 

The renormalization group treatment described in the 
next section can be set up within the original perturba- 
tion series (|43|) or the partially resummed series using the 
replacement ()46|) . Since we aim at a weak coupling ex- 
pansion, the partially resummed series makes only sense 
if the full Liouvillian is expanded in the same parame- 
ter as the renormalized vertices. Therefore, we will make 
use of the resummed series only at the end of the RG 
flow where perturbation theory in the renormalized cou- 
plings at a fixed physical cutoff scale can be used. For 
other problems like quantum dots in the charge fluctu- 
ation regime or systems in the strong coupling regime, 
it might be helpful to use the resummed series from the 
very beginning. 

Finally, we note some useful symmetry properties for 
the vertices and the Liouvillian (see Ref . \VA for the proof), 



G12 — — G21 , 
^12 ~ ^21 1 



TisLy^z) = 
Tr5Gi2 - 



reffi 



(57) 
(58) 
(59) 
(60) 
(61) 
(62) 
(63) 
(64) 



where 



Gil' = E^n' ' ^n' = E(^')ii' ' (65) 

p 



and the c transform A"^ of any dot operator A in Liouville 
space is defined by 



a: 



(66) 



Properties ((6T|) and (|62|) are important to show in time 
space that the reduced density matrix of the dot stays 
hermitian and the current stays real. Property ([5^ leads 
to conservation of probability, i.e., the normalization of 
the reduced density matrix stays constant. From this 
property it also follows that L'^J^ [z) has an eigenvector 
with zero eigenvalue: 



L'^^^{z)\x,{z)) = 



. (67) 

This eigenvector corresponds to the stationary state for 
z — > iO"*" and depends on the physical system under con- 
sideration. In contrast, the corresponding left eigenvector 
is unique and, according to Eq. ([M)) . is given by 



{xo{z)\ss') 



(68) 



As a consequence, in combination with property ([50]) . we 
obtain zero if the left eigenvector for zero eigenvalue acts 
from the left on the vertex averaged over the Keldysh 
indices: 



(2:0(2)1 G12 







(69) 



Therefore, by decomposing the vertex according to 

G\l = ^(Gn'+pGnO , G,v = Y^P^Tv , (70) 

p 

we see that the zero eigenvalue of Ll/-^{z) can only oc- 
cur in the resolvents when the part Gn' of the vertex is 
standing right to the resolvent. Therefore, to avoid this 
zero eigenvalue in the RG treatment, we will first use a 
certain perturbative treatment to eliminate the part Gn' 
of the vertex from the very beginning. This is described 
in the next section. 

B. RG equations 

First RG step. The first discrete RG step consists in in- 
tegrating out the symmetric part ^[fa{(^) + fa{~^)] — \ 
of the Fermi function in the contraction (|33p . This part 
depends only weakly on the frequency and creates no 
logarithmic divergencies in perturbation theory. Further- 
more, as explained in detail in Ref. [itI. it is the symmetric 
part of the Fermi function which allows the zero eigen- 
value of the effective Liouvillian L'g^ (E) to occur in the 
resolvents between the vertices. This part should be in- 
tegrated out before starting the continuous RG in order 
to show that the renormalization of the vertices is cut 
off by relaxation and dephasing rates. To get rid of the 
symmetric part, one decomposes the contraction (j33p ac- 
cording to 



71 



pp 



7i 



7i = p{Co) 



(71) 

(72) 
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1 1' 



1' 1 



1 1' 



1' 1 



1 2 



2 1' 



FIG. 3: The lowest order diagrams for the effective Liouvil- 
lian and the effective vertex when the symmetric part of the 
contraction is integrated out. s (a) denotes the symmetric 
(antisymmetric) contraction 7^ (7")- 



u 

1 2 



ii li li u 



2 1 12 



2 3 



3 1 12 



2 1' 



FIG. 4: RG diagrams for the renormalization of the Liouvil- 
han in 0{G^) and 0{G^) (first two diagrams) and the one-loop 
renormalization of the vertex in 0{G^) (last diagram). The 
slash indicates the contraction where the Fermi function has 
to be replaced by —dA 



d/2 



with tD = rjoj. Using this decomposition in Eq. 
one finds that each diagram, which is irreducible with 
respect to the full contraction 7, decomposes into a se- 
ries of blocks which are irreducible with respect to the 
symmetric part 7* (i.e., any vertical line hits at least one 
symmetric contraction) and connected to each other by 
antisymmetric contractions 7". The blocks which are ir- 
reducible with respect to 7** can be formally resummed 
into an effective Liouvillian L°-{E) or into effective ver- 
tices G°'{E), which obtain an additional energy variable 
to account for the reservoir contractions which cross over 
the effective quantities. The lowest-order diagrams for 
and are shown in Fig. [3l The first two diagrams 
correspond to the effective Liouvillian (no free lines) and 
the third one to the effective vertex (since two lines are 
free). Using the diagrammatic rules (1^5]) together with 
Eq. (|7ip and the convention ([36)) . we obtain for the first 
two diagrams. 



1 



L 



(0) 



and for the third one (including the interchange 1 ^ 1') 
, . „„„ 1 



P 72 ^12 



E12 



L 



(0) 



Gil, - (1 - 1') 



We use here the original perturbation series so that 
the unperturbed Liouvillian Lg-* occurs in the resol- 
vents. Performing the frequency integrations and omit- 
ting terms of order 0{1/D), we obtain the following per- 
turbative result for the effective Liouvillian and the ef- 
fective vertex containing the symmetric part of the con- 
traction: 

L^(i?) = Ls + Y.'^iE), (73) 

2 

-^^G,v{E,r - 4°))Gi,i , (74) 

Giy = Gil' ~ * (^Gi2 G21' — Gi'2 G21 j .(75) 

Analog equations hold for the effective current kernel 
S°(£') and for the effective current vertex 1^'°-. These 
are obtained by replacing the first vertex G by the cur- 
rent vertex r in Eqs. JTH) and ((75|) . 



1 Jl Ir^l l^i 11 



2 3 11' 3 2 12 1'3 3 2 23 3 1 1'2 

FIG. 5: RG diagrams for the renormalization of the vertex in 
O(G^) (two-loop). 



After integrating out the symmetric part of the Fermi 
function in this way, we obtain a new diagrammatic se- 
ries for the kernels analog to Eq. but the Liouvillian 
and the vertices have to be replaced by the effective ones 
and the contractions between the effective vertices con- 
tain only the antisymmetric part 7°. Furthermore, since 
the effective quantities have become energy dependent 
(also the effective vertex G° becomes energy dependent 
in higher order perturbation theory), one has to replace 



1 



E + Xi-L 



(0) 



G 



1 



E + X, ^ LI-iE + X,) 



G'^iE + X, 



in Eq. (|43p . Since the antisymmetric part of the con- 
traction ([7T|l does not depend on the Keldysh indices, 
only the effective vertex G" averaged over the Keldysh 
indices occurs in the new perturbative series. As a conse- 
quence [see Eq. (f69| ]. the zero eigenvalue of the effective 
Liouvillian can no longer occur in the denominator of the 
resolvents. 

Second RG step. The task of the second continuous 
RG procedure is to integrate out the antisymmetric part 
of the Fermi distribution function step by step. In each 
infinitesimal step, a small energy shell is integrated out 
and is incorporated into renormalizations of the vertices 
and the Liouvillian. However, instead of integrating out 
the energies on the real axis, it has turned out to be 
more efficient to integrate out the Matsubara poles of the 
Fermi distribution function on the imaginary axisji^i^^ 
i.e., in each RG step one integrates out one Matsubara 
pole starting from high energies. To obtain a contin- 
uum version at finite temperatures, one introduces a for- 
mal cutoff dependence into the antisymmetric part of the 
Fermi distribution by 



1 



.(A-KI) 



(76) 



where = (2n + l)TTTa are the Matsubara frequencies 
corresponding to the temperature of reservoir a, and 



2-kT 



for \uj\ > ttT 
for L < ttT 



(77) 
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is a theta function smeared by temperature. For A = oo, 
Eq. ([7S|) yields the full antisymmetric part /^(w) — 5 of 
the Fermi distribution. In each RG step, one reduces the 
cutoff A by dA, and integrates out the infinitesimal part 

ft - ft~'^^ = of the Fermi distribution. The new 

effective Liouvillian and the new effective vertices at scale 
A - dA, 

Lt'\E) ^ LUE)~dL'^s{E), (78) 
G'l~,'\E) = GUE)-dG^AE), (79) 

can be calculated technically in the same way as for the 
first discrete RG step. The only difference is that an in- 
finitesimal small part is integrated out so that the RG 
diagrams contain only one contraction involving the part 



dh.-^. Furthermore, since the diagrams should be irre- 
ducible with respect to this part, this contraction must 
connect the first with the last vertex of the diagram. 

Up to 0{G^) (which we call two-loop here^), the RG 
diagrams for the Liouvillian and the vertices are shown 
in Figs. [3] and [5l Using the definition 

li = P{^)ft{^) , (80) 
together with the convention Ls = L5, Gw = G^y and 



El...n + — Ls{Ei,,,n + UJl...n) 

we obtain the following RG equations: 



dLsjE) 
dA 



drfl 
dA 
drfl 
dA 



+ W12) 

I2 I3 Gi2{E) G2z{Ei2 + CJ12) G-3,l{Ei3 + UJis) 



(82) 



for the Liouvillian, and 
dGir{E) 



dA 



-^Gi2{E)Ui2G2r{Ei2+u;i2) - (1-1') 



I3 G23{E) 1123 Gil' {E23 + 1^23) nii'23 G32(-Ell'23 + '^11'23) 

f^^A ^ _ _ _ ^ 

73 G'i2(£') ni2 Gi'3(i?i2 -I- <I'12) nii/23 G'32(i?ll'23 + Wll'23) — (1 — 1) 

^73^G23(£;)n23G3i(£;23+W23)ni2Gi,2(£;i2+<:^12) - (l-l') 



(83) 



r 



for the vertex. Similiar RG equations hold for the current for the resolvents and the vertices, where Lg{E) and 

G^[E) are the renormalized quantities at scale A [in 
higher order in the couphng, also vertices G^ „(£■) with 
more than two indices can be generated]. If the exact 
RG equations in all orders are used, this perturbative se- 
ries gives the full kernels at each scale A. Therefore, it is 
possible at each step of the RG to resum all closed subdi- 
agrams between two vertices, leading to the replacement 



kernel Ti^{E) and the current vertex IJy{E) by replac- 
ing the first vertex in all terms on the r.h.s. of the RG 
equation by the current vertex. The initial conditions of 
the RG equations are given by Eqs. ([75|) - ([75)) from the 
first discrete RG step. Since ■Ji^^ — 0, the final solution 
at A = provides the result for the effective Liouvillian 
and the current kernel 



Ll"iE) - LsiE)U^o 



(84) 
(85) 



from which the reduced density matrix and the current 
can be calculated in Laplace space via Eqs. ([55]) and 

We note that one can stop at each step of the RG and 
use the perturbative series with the contractions 
at scale A together with the replacement 



1 



E + X, 



-(0) 



G 



1 



E + X, 



L^E^ 



G^iE + X,) 



1 



1 



E 



X,-L^{E- 



X,, 



E + X,- Ll/^ (E + X,) 



(86) 

where L'g'-^ (E) is the full effective Liouvillian at the end 
of the RG flow at scale A = 0. Using the perturbative 
series in a certain approximation, one can set up a self- 
consistent equation for Ll/^{E) at each step of the RG. 
However, this is only possible if the perturbation theory 
in the renormalized coupling is well defined. We will 
see that this is only possible at a certain scale A^, where 
some physical cutoff scale is reached, see Sec. Ill CI Up to 



11 



this scale, we will always use the renormalized Liouvillian 
Lg{E) in the denominator of the resolvents. 

RG in Matsubara space. Using the fact that the re- 
solvents and the vertices on the r.h.s. of the RG equa- 
tions are analytic functions in all frequencies Cui in the 
upper half of the complex plane, all frequency integra- 
tions can be calculated analytically by closing the con- 
tour in the upper half of the complex plane. The only 
poles occurring there are the poles of the contractions 
and their derivatives, given by 

7^ = -p(c^)T„^-^-i— 0tJA-K|),(87) 

uA ZTT \uj — «At„ + * At„ / 

denotes the Matsubara frequency cj" which lies clos- 
est to the cutoff A. After performing the integration 
we find that, due to the presence of the cutoff function 
P{lD) = jj-P^^i , the r.h.s. of the RG equations gives a 
negligible contribution for A 3> -D. Therefore, we can 
start the RG at Aq ~ I? and omit the cutoff function 



p((D) (finally, the precise ratio between Aq and D is de- 
termined such that no linear terms in D are generated, 
see below). As a consequence, only the Matsubara poles 
of the Fermi function in the upper half of the complex 
plane will contribute to the frequency integrations and 
all real frequencies are replaced by Matsubara frequen- 
cies. From now on, we write the frequency dependence 
explicitly and define the analytic continuation of the Li- 
ouvillian and the vertices in imaginary frequency space 

by 

LsiE,uj) = LsiE + iuj) , (89) 

^>)lfl:')2l^2(-^:^:'^l;^2) — Gl2{E + iU!)\ ,(90) 

where = lj", = w"^' correspond to Matsubara fre- 
quencies. With the definition 

n^^'^) = IF^- — \ (j, ^ ' (91) 

E + IUU — Ls(E,uj) 
the RG equations in Matsubara space can be written as 



^^^dA ' ^ G'i2(£',Ci^,AT„^,t^2)n(£'i2,AT„^ + + ^2) G2i(£'i2, At^^ + + ^^2, -^2, -At„ J 
-iGi2Tl{Ei2,AT^^ +UJ + uj2)G23TliEi3,AT^^ +uj + uj3)G^i , 



dGii'{E,u;,u;i,uj[) 



i {Gi2(£',w,wi,AT„jn(£'i2, At„2 + uj + uji)G2v{Ei2,At^^ + uj + uji,-At^^,u)[) 

-|-G23n(£'23, At„2 + UJ + UJ3) Gil' Il{Eii'23, -^T^^ + UJ + UJl + Uj'l +UJ3)G32 

- {Gl2Il{El2,■^-T^^ + u + ui) Gi>3ll{Eii'23, -^T^.^ + UJ + ui + uj[ +uj3)G32 - (1 ^ 1')} 

+ {G23Ti{E23AT^^+^ + (^3)G-3iU{Ei2,AT^^+0J + UJi)Gr-2 - (1^1')} ■ 

I 



(92) 

(1-1')} 

(93) 



In these equations, the index 1 = rjiy includes no longer 
the frequency variable, and we implicitly sum over all in- 
dices and Matsubara frequencies on the r.h.s. of the RG 
equations which do not occur on the left hand side (l.h.s.). 
Only positive Matsubara frequencies smaller than the 
cutoff A are allowed and each sum has to be written as 



2TTTa Ot^ (A 



(94) 



which reduces to an integral Jg duj for zero temperature. 
The frequency arguments of the vertices in the terms of 
O(G^) in Eqs. ^ and ^ have been omitted since 
they are not needed for the weak coupling analysis up to 
two-loop order, see below. 

The RG equations in Matsubara space are the final 
result of this section and are the starting point for the 
analytical solution in the weak coupling regime presented 
in the next section. Similiar RG equations hold for the 
current kernel and the current vertex in Matsubara space 



by replacing the first vertex in all terms on the r.h.s. by 
the current vertex. Using Eqs. and the 

effective Liouvillian and the current kernel follow from 



E^{E) 



Ls{E,oj 
i:~f{E,uj 



0)|a=o 
0)|a=o 



(95) 
(96) 



Finally, we note that all symmetry properties stated 
in Eqs. ([ST jl -lfM ]) are preserved under the RG flow (see 
Ref. (13 for the proof). 



Gl2{E,U!,UJl,UJ2) 
I^^(E,UJ,UJl,UJ2) 

Trs Ls{E,uj) 
Tis Gi2{E,uj,uji,uj2) 
Ls{E,u;r 
^^{E,ujr 

Gi2{E,UJ,UJi,UJ2y 
i^2{E,UJ,UJl,UJ2y 



- G2liE,LU,UJ2,UJl) 

- l]l{E,UJ,UJ2,UJl) , 



0, 
0, 



^Ls{~E,Lu), 
'^i{~E,uj) , 

- G2-l{-E,LO,L02,l^l) 



(97) 
(98) 
(99) 
(100) 
(101) 
(102) 
(103) 



■%{-E,L0,L02,UJl), (104) 
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where all energy variables are real. 

Similiar to the discussion at the end of Sec. Ill Ai prop- 
erties (Uni) and (jlOOp have the consequence that we obtain 
zero if the left eigenvector {xq{E,uj)\ of the effective Li- 
ouvillian Lg{E,ijj) for zero eigenvalue acts from the left 
on the vertex 



with 



{xoiE,co)\Gu = 



(105) 



compare with Eq. (l69|) . Therefore, it is not allowed that 
the zero eigenvalue of Lg{E,uj) occurs in the resolvents 
on the r.h.s. of the RG equations, proving that the RG 
is always cut off by relaxation and dephasing rates (this 
property holds in all orders of perturbation theory and 
within all truncation schemes). As we will see in the 
next section, this property is essential to prove that a 
systematic weak coupling analysis can be carried out in 
the generic case. 



C. Two- loop analysis 

In this section, we solve the two-loop RG equations 
and ([55)1 analytically in the weak coupling regime up 
to two-loop order. Weak coupling is defined by the condi- 
tion that the rcnormalized vertices Gi2{E, to, ui, UJ2) stay 
small compared to one throughout the RG flow so that a 
systematic expansion is possible on the r.h.s. of the RG 
equations. This condition is fulfilled if the various cutoff 
scales occurring in the resolvents n(£', uj) are larger than 
the energy scale Tk at which the vertices would diverge 
in the absence of any cutoff scales (the so-called Kondo 
temperature for the Kondo model). From the form of the 
RG equations, we see that the resolvents at scale A have 
the form 



with 



El.. 



(106) 



(107) 



and positive Matsubara frequencies < oj^ < K for 
/c = 2, . . . , 71. Here, E-\-iu \s the original Laplace variable 
at which wc want to calculate the final effective Liouvil- 
lian Llf-^ {E + iuj), and the index 1 corresponds to the 
contraction connecting the first with the last vertex of 
the RG diagram (where the Matsubara frequency is re- 
placed by the cutoff At^ ) . 

Expanding the resolvent (|106[) around its poles analog 
to the discussion at the end of Sec. Ill Al we arrive at the 
approximation ([M]) which contains the most important 
terms leading to logarithmic enhancements 



n^z) 



|xf(zf))(xf(zf)|, 



(109) 



[UJ + UJ2 + ■ 



where — ~ iF^, F^ > 0, denote the positions of 
the nonzero poles of the resolvent (we assume single poles 
here, but the following discussion holds also for other 
cases; note that the zero pole of the stationary solution 
can not occur in the resolvent as discussed at the end 
of the last section) . Since all Matsubara frequencies and 
the relaxation/dephasing rates are positive, we see that 
the resolvents can not become large. Using < Wfe < A 
for k = 2, . . . ,n, and setting cj = 0, we find that the 
resolvent is cut off at the scale 



A ~ max{T,,|Si...„-/if|,ff} 



(110) 



where the maximum is taken over all values of the oc- 
curring indices. Here, temperature is a trivial cutoff pa- 
rameter, because, for A < 27rTa, the sum ([94)1 over the 
Matsubara frequencies for reservoir a reduces to one term 
n = and the cutoff Aj-^ = ttTq becomes independent 
of A. Therefore, temperature is a unique cutoff for all 
terms on the r.h.s. of the RG equations, like in equilib- 
rium problems. This trivial cutoff is set to zero in the 
following, i.e., = 0, and we discuss only the nontriv- 
ial dependence on the other cutoff scales. The minimal 
cutoff scale occurs for Ei...n — hf = and is given by 
the relaxation or dephasing rates Ff . These points define 
the positions of resonances where renormalization-group- 
induced logarithmic enhancements or suppressions have 
to be expected. However, as we will show in the fol- 
lowing, these logarithmic terms can be calculated sys- 
tematically by perturbation theory in the rcnormalized 
couplings, provided that the weak coupling condition 



Ac EE max{|£:|, > Tk 



(111) 



(108) 



is fulfilled, where by convention hi = h^~^ denotes the fi- 
nal (physical) rcnormalized oscillation frequency at scale 
A = (we will see that the difference between hf and 
hf^^ is proportional to the final rcnormalized coupling, 
i.e., only a small perturbative correction). Ac is an im- 
portant energy scale separating two energy regions where 
the RG equations are solved in a different way. It is given 
by the maximum of the Laplace variable E, the chemi- 
cal potentials of the reseroirs (giving some voltage V), 
and the oscillation frequencies hi of the different phys- 
ical modes (e.g., the rcnormalized magnetic field in the 
Kondo problem) . This is roughly the maximum value the 
various cutoff scales — h^\ = \E + J2k Vk^^a^, — 

of the resolvents can take, see Eq. (|110p . Thus, for 
A > Ac ^ Tk, the cutoff scales do not play an important 
role and we get Ja ^ 1, where Ja is the order of mag- 
nitude of the vertex at scale A. Since A is the relevant 
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energy scale in this regime, the order of magnitude of the 
relaxation/dephasing rates is given by ~ AJ| <C A 
(note that the RG for the LiouvilUan starts in second 
order in J). Therefore, the relaxation and dephasing 
rates are small perturbative corrections in the denomi- 
nators of the resolvents and do not lead to any cutoff 
of the RG in the regime A > Ac. Since all vertices are 
small, we can systematically truncate the hierarchy of RG 
equations and expand the solution systematically around 
the poor man scaling solution (i.e., the lowest-order so- 
lution for the vertex in the absence of any cutoff scales 
Ac = Tq, = 0). This gives a certain initial condition 
for the RG at A = Ac presented as a power series in 
Jc = Ja=a^. 

In the second regime < A < Ac, the RG for Tf is very 
weak and will be roughly cut off by Ac. The reason is that 
there are many terms on the r.h.s. of the RG equation 
involving different cutoff scales |£'i...„ — but usually 
one of them will be given by Ac already in second order 
in J. In this case (see a comment below for the other 
cases), we get Tf ~ AcJc for all < A < Ac because 
becomes smaller for decreasing A. This means that 
even at resonance, the minimal cutoff scale is given by 
r ~ AcJ^. Therefore, by expanding the solution for the 
vertices systematically in Jc, we get in the worst case at 
resonance a series of the schematic form 



Ja ^ Jc { i + Jc In 



Ac 



\A + iT\ 



(112) 



For A ^ 0, the logarithmic term becomes maximally of 
the order ~ Jc In Jc ^ 1, which is a perturbative correc- 
tion in the weak coupling case Jc <Si 1- Therefore, un- 
der condition (|llip . we stay in the weak-coupling regime 
and the RG equations can be solved perturbatively in 
Jc in the whole regime < A < Ac. We note that this 
fact relies essentially on the condition that all resolvents 
on the r.h.s. of the RG equations contain some relax- 
ation/dephasing rate Tf. As explained at the end of the 
last section, our RG approach gives this property in the 
generic case in all orders of J by construction. 

In case that the RG equation for contains a smaller 
cutoff scale Ac <C Ac in second order in J, we expect 
~ A^(J^)2 for A < A^, with J^ = Ja=a^. In this case, 
the logarithmic term in Eq. (|112p leads to contributions 

Jc In for A ^ 0, giving rise to additional enhance- 
ments and sharper features at resonance. However, even 
for AJ, — + 0, there is no divergence since the cutoff scale 
Ac will certainly occur in some higher order term on the 
r.h.s. of the RG equation for F^. Thus, the minimal cut- 
off scale will be of order ~ Ac Jc with k > 2. This gives 
a maximal value ~ fc Jc In Jc for the logarithm which is 
again a perturbative correction just enhanced by a factor 
of k. This shows that the height of logarithmic enhance- 
ments at resonance are expected to be increasable only 
by factors of 0(1), but the sharpness of features at reso- 
nance (which are controlled by F^) can become orders of 
magnitude smaller. 



As a consequence, we have seen that for < A < Ac, 
we can perform a perturbation theory in Jc, which is 
the order of the vertex at scale Ac. This means that we 
can equivalently stop the RG at A = Ac and use the 
perturbative series ([43]) with the contraction 7^^^ at scale 
Ac together with the replacement 



1 



E + Xi- L 



(0) 



G 



E + X, 



L^^{E + X, 



for the resolvents and the vertices. In contrast to the 
perturbative series at scales A :§> A^ this perturbative 
series at scale Ac is well defined and can be used alter- 
natively to the RG approach. Furthermore, as explained 
in Sec. Ill A[ the perturbation series can be partially re- 
summed, leading to the replacement 



1 



1 



E + X,, 



Xi) E + X^~ Ls"{E + X^) 



(113) 

i.e., the final full effective Liouvillian can be written in the 
denominator. This series has the advantage that the os- 
cillation frequencies hi, defining the resonance positions 



El. 



E 



k 



(114) 



and the relaxation/dephasing rates Fj, cutting off the log- 
arithmic enhancements at resonance, are the final phys- 
ical ones at scale A = 0. This is expected on physical 
grounds and leaves no question open what the precise 
prefactor of these energy scales is. Using the replace- 
ment (jll3p . one can cither write down directly the per- 
turbative series or one can use the RG equations [again 
using the replacement (|113p to define the resolvents] and 
solve them perturbatively in Jc. Both approaches give 
the same because the RG equations are formally exact. 

Having shown that a weak coupling analysis is well 
defined for all cutoff scales under the condition (|llip . we 
proceed to show analytically the perturbative solution of 
the RG equations in all details for the two regimes A > Ac 
and < A < Ac in the generic case. 



1. RG above Ac 

Lowest order. For A > Ac, we define the reference 
solution for the vertex by considering only the first term 
on the r.h.s. of RG equation ([93]) with H(i?i2, A uj + 
uoi) i.e., by setting all frequencies to zero and 

omitting the cutoff scales from E, /ii2, and Lg{E). This 
defines the leading order RG equations 



dG^, 
dA 


1 
A 




Ml) 


-(1^1' 


"-'ii' 
dA 


1 
A 




' Ml) 

^21' 


- (1 4-> 1 



■)} 



(115) 
(116) 
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The initial condition for these RG equations is the bare 
vertex. The order of magnitude of the leading order solu- 
tion is denoted by the dimensionless parameter J ~ • 
The connection to conventional poor man scaling is es- 
tablished by recognizing that the leading order vertices 
have the same form in Liouville space as the original ver- 
tices, given by Eqs. ([TU)) . (IHU]) . and (PTjl . i.e., one can 
prove thatjii 



Ml) 


= [512, ■]- , 


(117) 


Ml) 


= [512,-]+ , 


(118) 


f7(l) 
-'12 


— ''12 '-^12 I 


(119) 



with 



c^.9ii' 1 r I. 



1')} 



(120) 



Thereby, the form (jll9p of the current vertex can only 
be proven if one takes the Tr5 over the local quantum 
system, i.e., it holds only for the combination Tvsl22^ ■ 
Implicitly, for all following equations, we will always con- 
sider this combination for the current vertex because this 
is finally needed for the calculation of the average of the 
current. Equation (|120p is the usual poor man scaling 
equation which can also be derived on a pure Hamilto- 
nian level by, e.g., leaving the t matrix invariant.^ 

Next we set up the lowest order RG equation for the 
Liouvillian Ls{E, lo) by considering the first term on the 
r.h.s. of Eq. ([5^ and replacing the vertices by the leading 
order ones. Furthermore, we replace Ls{E,uj) by L^g^ in 
the resolvent. This gives 



dLs{E,Lo) 



dA 



L 



(0)^ 



Ml) 

'^21 



with 



/Ca(z) 



In 



/2A 



\ A-iz 



(121) 



(122) 



To extract the lowest-order term ^ J from this equa- 
tion we treat the terms ^ z/ A of /Ca(z) separately by the 
decomposition 



iz 
2A 



/Ca(z) = /Ca(z) -f 
so that JC\{z) is integrated by the function 
JCa{z) 



Fk{z) 



A In 



2A -iz 



A - iz 
iz / (2A-iz)A 
~ 2" \ 2{A -izf 

with the following asymptotic behavior: 

, z , 



Fk{z) 



A 



ln2 + 0(^)2 



for A > \z\ 



(123) 



(124) 



(125) 



(126) 



Using Eq. (|12ip . the second term on the r.h.s. of 
Eq. (|123p leads to the following RG equation for the Li- 
ouvillian in leading order: 



dL'^g\E,Lj) _ 1 



G[WEi2+iu;-L'^')G. 



dA 2A 
with the initial condition 
(1)^ 



(0)x Ml) 



L's'iE,u;)\A=Ao = 



(127) 



(128) 



When integrated, we obtain a contribution L^g'' (E, lo) ~ 
J, i.e., one power less than expected due to the factor 
[compare with RG equation (|115p . where the same 
happens]. In contrast, the contributions from /Ca(-z) are 
of second and third order in J, as will be discussed below. 
We write 



L«(£;,c.) 



's 



with L 



(1) 



L'-g' - {E + ioj) (129) 
L^s\f = 0,w = 0), and 



dL 



(1) 



dA 

dZW 

dA 



^GW(M12 



-(o)-.Mi) 



=.(1) 



Ml) 
'^21 



(130) 
(131) 



We note that ^^^(i) can be interpreted as the Z factor 
in Liouville space at scale A, and it can be shown that 
L^g"^ and Z^^^ are hermitian operators. Similiar equations 

can be set up for the current kernel E^"'^^ (iJ, to) in lead- 
ing order by replacing the first vertex on the r.h.s. of 
Eqs. (|130p and (|13ip by the current vertex. 

When integrating the RG equations (|130p and p3ip up 
to Ac, we obtain a linear contribution in the renormalized 
coupling Jc at scale Ac. Since the perturbative treatment 
in Jc for the regime < A < Ac can only give corrections 
~ Jc to the Liouvillian, we know that the final effective 
Liouvillian up to 0(Jc) is given by 



iE,L0) 



(132) 



.'i^' - {E + tiu)Z^'^'' 
and Z(i)^ = Z(i)| 



O(Jc') 



with L^g''" = L'g'' \a=a^ ana ^---'^ = zi^-'|A=A, 

Second order. With the vertex and the Liouvillian in 
leading order, we can now expand the full RG equations 
systematically around these reference solutions and cal- 
culate the higher orders. The Liouvillian and the vertex 
are written as an expansion in J, 



LsiE,iu) 



- r(0) 



+ 4'\e,cu) + L],^'{E,u;) 



Gi2{E,UJ,UJi,UJ2) 



G 



(1) 

12 



+ G'f^{E,LO,iO^,L02) + 



(133) 



(134) 



where h'g'^ , g'^^ 



J with k < n, see below). 



J" (possibly with additional factors 
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To get an RG equation for L^g\E,uj) and 

Gi2 {E^bj,uJi,uj2), we insert the expansions (|133p 
and (fm]) into the r.h.s. of the fuU RG equations ([92]) 
and For the RG of the LiouviUian (vertex), the 

resolvents are expanded such that we collect all terms 
in 0[{^fj^] {0[l(f with k > 1, and 0{^J^) 

[0{j^J^)] on the r.h.s. of the RG equation, where A 
is some cutoff scale arising from £12, Ls{E), or the 
frequencies. To achieve this, we leave out L^g\E,uj) in 
the resolvents [leading to terms of 0{J'^) on the r.h.s. of 
the RG equation], and use Eq. (|129p . This gives for the 
resolvent 

1 

iI{{E,uj) « 



iE + iL 



(0) 

s 

1 



iL 



(1) 



+ (a; - iE)Z^^) 
1 



z(i) 



uj — iE + iL 



(0) 



iL 



(1) 



1 



2 



(135) 



with the hermitian operator 



/(I) 



1 



(1) 



(1) 



Z(i) 



2 



4"^ z(i)) 



(136) 



equations). For the resolvent in the first term on the r.h.s. 
of this RG equation, we expand (|135p in the following 
way: 



zn(Si2,A 
i-z(i 




(137) 



Together with the two vertices, this approximation con- 
tains systematically all terms of 0{j^J^), 0{j^J^), and 
0[^(-^)'^ J^], with fc > 1, which are important to calcu- 
late the vertex up to 0{J^). The last three terms on the 
r.h.s. of Eq. ([QS]) are already ~ J^, therefore we replace 
the resolvents by their lowest order term ^ j^: 



iIl{E23,^ + UJ + iU3) 

i n(i?ii/23, A + UJ + uji + u;[ + LJs) 
in{Ei2,A + uj + uji) 



1 



A + LU3 

1 

A + CJ3 
1 



A 



Vertex in second order. We start with the RG equation 
(j93p for the vertex (or equivalently the current vertex 
by replacing the first vertex by the current vertex in all 



Note that uj^ is an integration variable and has to be 
kept in the resolvent. As a result, we get the following 
RG equation for the vertex in second order in J: 



dGfi!(£:,w,wi,wi) 



duj3 



dA 



G 



(1) 
12 



A 



d , A + UJ + ijji — iEi2 + iL a 

-jrln 7 ^ 

dA A 



(0)- 



<-^21' 



1') 



\{G''^Gf^,{E^2,A + uj + uj,,-A,uj[)^G'^^^{E,uj,uj^,A)G\^^, - (1 ^ l')} 



1 



A + UJ3 



Ml) r<(i) f=<('i-) 

■-^23 '-^ll'^32 



— / dcj3 



1 iMi)r{i)r(i) 

— 1^12 (-l'3^32 



A + UJ3 



G^G^G^-il^V)} 



(138) 
(139) 
(140) 



r 



The first term [Eq. (|138p ] on the r.h.s. containing the 
logarithm induces the frequency dependence of the ver- 
tex. If we neglect in this term the derivative of G^^-* with 
respect to A (giving rise to terms of 0[j^{^Y J^], with 
/c > 1, contributing to G'^'), we can decompose the ver- 
tex in the following way to solve the above RG equation 
in 0(J2): 



with 



G^^^){E,uj,uj^,Lj[) 



(0) 



^(1), A + uj + uji- iEi2 + iL'g - ^(1) 

<^12 in - Grjl' 



In ^^±^i+^A_^El!l±lIllG^^) . (142) 



Using this solution in the second term [Eq. (I139p ] on the 

G^^y{E,Lj,LJi,Lj[) = G^j^^'' + G'"^-^\e,uj,uji,uj[) , (141) r.h.s., and neglecting again terms of 0[-^(-^)'^J^], we can 



i(2a) 



^(26) 
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use the approximations in Eq. (|139p . which lead to two terms canceUing precisely 

the second term of Eq. (|140p . Thus, we obtain the fol- 
G^y{Ei2^k + ijj + uJi,—k,tjj'i) w —\n2G^y\G'i^^ lowing RG equation for the frequency-independent part 

-(2a) 



G(f (£;,^,^i,A) « -ln2GWG(!) ^n' of the vertex in 0( J 



1^ = i{G«Gr;^+Gir)G«-(i^i')} 

- i {g«z(^)G« - (1 ^ 1')} - ^G«G«GW , (143) 



where the initial condition is given by the second term of 
the inital condition (|75|) for the vertex. It can be shown^-^ 
that the form of the initial condition is preserved [with 
the vertices given by the leading-order solutions (|117p 
and (|118p ] if one considers only the first term on the r.h.s. 
of Eq. (|143p . Therefore, we decompose the frequency- 
independent part of the vertex in second order as 

Gpi?' = iGfi?^' + Gfi?^) , (144) 

with 

^(2ai) _ /-(I) Ml) /-...x 

and G^^?^'' fulfils the same RG equation (|143p as G'f^\ 
but with zero initial condition. Since no explicit imagi- 
nary factors occur in Eq. (I143P , the decomposition p44p 
can also be viewed as a decomposition of the vertex into 
real and complex parts (provided there are no complex 
terms in the initial condition for the original quantities). 
Therefore, the two parts have a completely different phys- 
ical meaning. Whereas the part G^^?^'* is even impor- 
tant to calculate the rates in second order in J (see be- 
low), the part G^^?^'* denotes a renormalization of the 
coupling constants in two-loop order, which can lead to 
logarithmic corrections of the form ^ In , where Jo 
denotes the original coupling constant (see Sec. IIII B H 
where such terms are explicitly calculated for the Kondo 
model). Such terms arc not well defined in the scaling 
limit Jq ^ and, therefore, should be taken together 
with the leading order vertex by redefining certain char- 
acteristic low-energy scales (like the Kondo temperature 



for the Kondo model). Thus, in the following we will 
redefine the leading-order vertex by the replacement 

G^ii' G^iv + G^iv"^ : (146) 

and will consider only the part G^^?^'' explicitly in all 
equations. We note that it can not generically be shown 
that the replacement ()146p can be accounted for by just 
renormalizing the Kondo temperature in the lowest-order 
vertices. For the Kondo model, this can be shown in all 
orders^'^ because only one coupling constant remains in 
the scaling limit Jq — > 0, D — > oo, such that Tk stays 
constant. However, for more complicated models includ- 
ing orbital degrees of freedoms, interference effects, etc., 
it can happen that the matrix structure of G^^^?^^ in Liou- 

ville space is different from that of G[^! . In this case, new 
terms can finally arise from the second-order part of the 
vertices, and their influence might be quite nontrivial. 

We note that all equations also hold for the current 
vertex Jt(^) by replacing the first vertex in all terms on 
the r.h.s. of the RG equations by the current vertex. 
Thereby Eq. ()145p is only valid if the trace Trs over the 
local quantum system is taken from the left (which we 
always implicitly assume for the current vertex and the 
current kernel). 

Liouvillian in second order. We proceed with RG equa- 
tion ([^^ for the Liouvillian (or equivalently the current 
kernel by replacing the first vertex by the current ver- 
tex in all equations). Inserting the expansion (I134p into 
this RG equation, using Eqs. (fUTj) . (fTiS]) . and iMe]) . and 
neglecting all terms of 0{J^) on the r.h.s., we obtain 



dLs{E,uj) 



G[\^U{Ei2,A + uj + uj2)G'^^ 



(1) 



dA 



+ I G^l^ n(^i2 ,A + L0 + L02) G^f + I G\f" n(Si2 , A + CO + C^2) G^^' 



(2ai) - 



^(1) 



- G[\^ n(^i2, A + CO + C02) G'^{' (S12, A 

- UJ 



C02, -UJ2, -A) + G^f ^ {E, CO, A, C02) n(Si2, A + CO + C02) G^j 



(1) 



■iG[\^n{E,2,A- 



Lj2)G^i^U{Ei3,A + uj + uj3)G^^^ 



(147) 

(148) 
(149) 
(150) 
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For the terms (IMT]) . ((T48)) . and ((T50)) we need the Eqs. p35)) . (fT36l) . (fT22)) and (fT23)) . can be expanded as 
resolvent integrated over frequency, which, by using 



ij dL02TliEi2,A + ^+i^2) ~ h -— j /Ca(Si2+zc^-4°^-4'Mi-— j 
7(1) 7(1) 

« ^ (1 - — ) (£^12 + - 4°^ - ii'^) (1 - — ) + ^a(^i2 + - 4°^ - 4'^) 

« ^ + - ii"^) - ^ ((i?i2 + *c.)Z(i) + - Z(i)4) - Lfz^^^) + ICa{E,2 +iu:- 4^ (151) 

I 



The first term is of order O(^), the second one of or- 
der J), and the last one contains terms of 0(1) and 
0[(^)'^], with fc > 1. Therefore, when multiplied with 
J" and integrated, the first term gives a contribution of 
order 0(AJ"~^), the second one leads to 0(AJ"), and 



the last one gives 0(AoJ") or 0(AJ"). Thereby, the 
terms ~ Ag J" are cancelled by the initial condition from 
the first RG step (see below). Therefore, for the calcula- 
tion of -^L^g\E,Lj), the terms (IWl) and give rise 
to 



lA {'^i2^^A(i?i2 +^u^- Lf)G^^ } (152) 
+ ^ {g[\\e^2 + ^. - Lf)G^^^-^ + Gg^^\E,2 + - lf)G^ (153) 
- i {{E,2 + ^u)Z^^) + 4) - Z^^^lP - 4)Z(1)} G«. (154) 



As is shown in Appendix [Aj the other two terms, 
Eqs. (|149p and (|150p . have nearly no effect when ex- 
panded systematically, one just has to replace the func- 
tion Fa{z) in Eq. p52p by the function 



Fj,{z) = Fa{z) 



In 2 



(155) 



According to the three terms Eqs. (|152p - p54p . we de- 
compose 



Lf\E,c.) = L^f>iE,cu)+Lr +L 



- (2a 



(26) 



-(2c) 



with 



dL 



(2h) 



dA 



dZ^ 
dA 



dL 



(2c) 



dA 



dZ 



(2c 



dA 



-iG[\^Fj^{Ei2 + iuj 



Lf)G^^^, (157) 

rr<(l)/- ^(0)x^(2ai) 



-G 



(2ai) 



12 



L 



(0) 



(158) 



i r^(l)^(2ai) 

-2A 1^12 



^(2ai)^^± 

"-^12 ^oT 



51^} . (159) 



2A 



(1) 

'12 



{(/2l2^W 



4^- 



-zW4)-x, 



^G^zm^. 



G 



(1) 



(160) 
(161) 



Note that, according to Eq. (|126p . the initial value of 
L^g'^\E,LL)) is given by 

4")(i?,c.)|A=A„ = -zln(2)AoG^^2)G^y|A=A„, (162) 



{E + iuj) (Z'-^''^ + Z'^'^^^) , (156) which coincides with the first term of the inital condition 
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((74)) from the first RG step if we choose 



Ao 



161n(2) 



D. 



(163) 



One can group the various terms occurring in 
Eqs. (|156p ~ (|16ip regarding the occurence of the complex 
factor i. The RG equations (|160p and (|16ip do not con- 
tain an expHcit factor i and integrating them up to Ac, 
one generates a second-order contribution which is neg- 
ligible compared to the first-order terms ()130p and (I131|) 
(up to this value of the cutoff parameter, no logarithmic 
contributions are generated, i.e., the terms just have one 
factor more in J). In contrast, the second-order terms 
containing the factor i are not negligible because they do 
not occur in first order. They are generated by Eq. (jl57p 
[with F^{z) replaced by its real value] and by Eqs. (|158p 
and (|159p . This is a general rule that the terms contain- 
ing the factor i are always generated one order higher in 
J compared to those without this factor. 



2. RG below 



The two-loop RG until Ac has resummed all leading 
and subleading logarithmic contributions in In into 
the renormalized vertices. This means that we have con- 
sidered all terms of the form 



I'' In' — 
A. 



I 



1,/c 



(164) 



for the Liouvillian and the vertices, where Jq denotes the 
original coupling constant. The terms with I = k — 1 and 
I = k — 2 arc defined in our terminology as one-loop and 
two-loop terms, respectively, irrespective of the topology 
of the diagrams, which is not a unique definition and 
depends on the formalism used. Roughly speaking, at 
A = Ac, the band width D has been replaced by an effec- 
tive band width Ac and the bare coupling constant is re- 
placed by a renormalized one Jc (including one-loop and 
two- loop renormalizations). This eliminates all the log- 
arithmic contributions of Eq. p64p , and a simple power 
series in Jc remains, which is well defined for Jc ^ 1, as 
described in detail at the beginning of Sec. Ill CI 

As a consequence, we solve the RG equations pertur- 
batively in Jc in the regime < A < Ac. Furthermore, 
we replace the Liouvillian in the resolvents by the full 
effective Liouvillian L'^J^ {E^ui)^ and we use the approxi- 
mation l|54p 



RG equations, these corrections lead to terms of order 
0{J^) without any logarithmic enhancement, which we 
neglect in the following. The eigenvalues of Ls are given 
by hi — iTi with > 0. These eigenvalues are finally 
calculated self-consistently from Eq. ([5T|l . In contrast, 
the Z factor is expanded as 



Z 



1 + + 0(j2) , 



(166) 



where Z^^^ ^ Jc ^ 1. This first-order correction to the 
Z factor can in principle be calculated from the effec- 
tive Liouvillian up to first order in Jc, given by the re- 
sult (jl32p . If the eigenvalues of Lg'' are separated well 
compared to the first-order corrections to the effective 
Liouvillian, we get 



2(1) ^ _zWc 



(167) 



where zj^'''^ is the diagonal part of Z^^'^'^ with respect to 
the eigenbasis of However, we will see later that 

Z^^^ does not enter our final result so its precise value is 
of no relevance. 

In the RG equations, the variable E is replaced by 
Ei...n and oj is the imaginary part of the Laplace vari- 
able plus the integration variables (the cutoff A and some 
Masubara frequencies). Thus, in the end the low-energy 
cutoff will be given by an expression of the form 



iuj - Lc 



(168) 



where from now on uj is the imaginary part of the origi- 
nal Laplace variable. Resonance positions are defined by 
Ei...n = hi, see Eq. (|114p . At these points, logarithmic 
terms of the form 



J^ In' ^ 



; = fc-l,fc-2, 



,0 



(169) 



are generated with A = -Ei...„ — + iVi- At reso- 
nance, the logarithmic terms are of order In' Jc since 
Vi ^ ^cJc^ leading to enhanced, but still perturbative 
corrections. 

Our first aim is to collect all terms of the form (|169p 
with A: < 3 and I = fc — 1, fc — 2, i.e., all terms of the form 



0(1), 0(Jc), 0(J2), 0(J2 In 



Ac 



A-f zFl 



(170) 



z 



E + iuj ~ Ls 



(165) 



together with the convention ([55]) . The corrections to this 
approximation are at least of order 0( Jc) and contain no 
poles in the variable E + icu (at least if we assume that 
the projectors on the eigenvectors do not contain poles 
or do not contribute). Therefore, when inserted in the 



Only terms of 0(Jc) without a logarithmic factor are ne- 
glected [note that approximation (|165p for the resolvent 
is already neglecting such terms]. Therefore, it is not 
necessary to calculate the Liouvillian in third order in J 
for A > Ac since no logarithmic contributions are gener- 
ated up to Ac. Furthermore, for A < Ac, we see that the 
lowest order term on the r.h.s of the RG equation ([^^ for 
the Liouvillian is already of 0{J^). Therefore, we need 
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the vertices only up to O(J^), and it is not necessary to 
calculate the vertices up to third order for A > A^. 

Vertices. We start with the perturbative expansion for 
the vertices. Up to order 0{J^), we get from the RG 
equation ^ together with Eqs. pBS]) . p66)) . and pM]) 

dG {E,uj,uj I, u![) _ 
dA 

- {°"'° A + -('"'')}- <™) 

Integrating this equation from A = Ac to A, we obtain in 

oiJ!) 

where clll^iE ,uj,uji,uj'i) is the value of the vertex in sec- 
ond order at A = Ac. Inserting this value from Eqs. (|14ip . 

I 



(fT42l) . and (fT44)) . we find up to 0{ J^) 

+ {Ciy^ln A±i^^Lp^G«; - (1 ^ 1')} ,(172) 

where we have used that A12 = £'12 + iuj — £5 •* + 0{Jc)- 

Note that the part G^^^?^'' is always taken together with 

G^^-l, see Eq. (|146p . therefore it does not occur in 
Eq. dlZll). 

Liouvillian. We proceed with the perturbative expan- 
sion for the Liouvillian given by the RG equation (j92|) . 
In the second term on the r.h.s. of this RG equation, 
we can replace the vertex by the lowest-order term G^^^'^ 
and the resolvent by II(E,uj) — — - — according to 
Eq. (|165p . In the first term on the r.h.s., we have to 
consider the first-order correction Z'^^ of the Z factor as 
well, and we have to use the second-order term for the 
vertex given by Eq. (I172p . which gives 



G^^iE,.,A,.,) = ^Gg-^^^ + GW-lngA_!^G^i)- _ G^- In ^ + '^^^ G^c ^ (^73^ 
G^^{Ei,,A + ^ + ^,,-oj,,-A) = zG^f + Gg- In ~ G^- - G^f In ^ + G^- . (174) 

For the frequency integral over the resolvent, we use 

duj2—, T- = -z/Ca(A) = -i In ^5 ^ • (175) 

Using these replacements in Eq. ((92|) and collecting the various terms according to their order in Jc, we obtain up 
to 0(J3) 



with 



Using 



Ls{E,cj) = + 4^^(i?,c^) + L^s\e,u;) + l[^''\e, cj) + Lf\E,u;) , (176) 



'-^^^^ = -zGW^^a(A,2)gW% (177) 



dA 



-^Glf {z(i)^a(Ai2) + k;a(Ai2)z(i) } G«' 

+Glf /CA(Ai2)G(f + Gf/ ^)^/CA(Ai2)G(f , (178) 
^ -^G^^^^^,iAl2)G^^^ln'-^^^ 

-^G^^ r d., I ^-^G^f In ^^^^^ + In ^^^^^G^f ^--1 G«^ 

J A \^2~ lAi2 23 Ac Ac 23 - lAi3 J ^1 

+zGjf/CA(Ai2)G^f /CA(Ai3)G^f . (179) 



<^ In — — + In — — } = In — In — In — In — — 180) 

[x + a\ Ac J \ Ac J x + b Ac Ac Ac Ac 
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for the third term on the r.h.s. of Eq. (|179p . we obtain after some manipulations 

^ _2^^a)c L 2A-^A,.^(,. 2A_^A^ _ ^^^G^f In \ G^f . (181) 

dA \ Ac ^■^ Ac Ac 23 Ac J 31 

Equations (|177p and ()178|) can be integrated easily from A = Ac to A = by using 

/Ca(A) = Afa(A) , Fa(A) = FA(A) + ^ (^ln-|^ + l) , (182) 

where ^'^(A) is defined by Eqs. pTO]) and p^ . Using Fa=o(A) = -^ln2, we get 

dA/CA(A) = -F;JA) - ^ (^In^ + l) . (183) 

The value of the Liouvillian at A Ac in second order has to be taken from Eqs. (|156P " (|159p . When integrating 
Eq. (|177p . the part from the first term on the r.h.s. of Eq. p83p cancels term (I157P for L^g"''^'^ {E , lu) if we neglect the 
difference between Lg-* and Ls [leading to terms of order 0{J^)]. Thus, for A = 0, we obtain the contributions 

L«(£;,a;)A=o = 4'^'-(i? + »c.)Z(i)^ , (184) 

In— ^ + 1^4^' (185) 



4')(£;,a;)A=o = Lf^^ + Lf^^~(i? + zc.)(Z(2^)^ + Z(2^)^)-iG«^Ai2 (in^ 
Lf)(£;,^)A=o = -icif (zWAiJlnl^lUAi^flnl^lU^^^lcW^ 



12 

-ziclf Ai2 (inl^l) G^f )^-ziGf-'^Ai2 G^f • (186) 

Thereby, we have only included the logarithmic terms in L^g°'\ i.e., terms ~ without a logarithmic factor are 
neglected. 

Finally, integrating Eq. (|18ip . we obtain 

Lf\E,Lj)A=o = -*Ac dxG('2)^ln(x - i^) cg^' ln(x ~ ^^) G^j^= , (187) 

JO -^c Ac 

Neglecting again terms ~ J^, the double logarithmic integral can be replaced by (with a = —i^^ and b = —i^^) 

da; ln(a; + a) ln(a; + 6) w -ialn^ a- ifeln^ 6 + alna + 61n6+ ialn^ ^ - aln^ , for |a| < |6| , (188) 

Zi zoo 

and analog for \a\ > |6| by interchanging a ^ b. For some special cases, (|188p can be written as 

ax ln(x + a) in(x + t>) ^ <^ -alnaln6 for |a| « |6| - 1 ' 

Equations (|184p - (|187p together with Eqs. I|188p - (|189p are the final results for the effective Liouvillian (or the current 
kernel if the first vertex in all terms is replaced by the current vertex) for a generic model of a quantum dot in the 

Coulomb blockade regime. The vertices g[\]'' and Gf^?'''' at A = Ac follow from Eqs. ([TT^ and pi5)) . The first- 
order quantities ig^'^ and Z'^^^'^ follow from the solution of the RG equations (|130p and p3ip at A = Ac, and the 

second-order terms L^g''^'^, L^s'^^'^ , Z^^'')^, and Zf^^)^ are determined by the RG equations (|158p - (|16ip . 

Looking back at Eqs. (|177p - (|179p . we see that the third-order terms involving three vertices on the r.h.s. of the RG 
equations enter only explicitly via the last term on the r.h.s. of Eq. p79p . This part is of 0{J^) and does not contain 
any logarithmic contribution since -ftrA(Ai2) — > *^7J f'^'' A — > 0. This means that the logarithmic contributions are only 
generated by the terms in second order in the renormalized vertices, but including their corrections in second order in 
Jc via Eq. (|172p . Both the imaginary parts iG^2ai)c ^^^^ |.]^g frequency dependence generate logarithmic contributions. 
Implicitly, third-order terms in the vertices are also present in the last term on the r.h.s. of the two-loop RG equation 
p43p . which determines the contribution G^2i2)^ ggg (jl46l) . As already mentioned after Eq. (|146p . it may happen 



21 



for more complicated models than the Kondo model that this part can change the leading-order vertex in a nontrivial 
way (i.e., not only by a change of the Kondo temperature) leading to new logarithmic contributions not expected 
from the matrix structure of the leading order vertices. 

For u; = and away from resonance IA12I, IA13I ^ F, we have a — —ia and b = —ib, with a and b being real. In 
this case, we get for Eqs. (|188l) and (|189|) 

— if dx ln{x ~ ia) \n{x — ib) fa —a\n^\a\-\ — 61n^ 161 — aln lal — 61n l&l aln^ 1^1 

Jo 2 2 2 6 

-i^f|a|ln|a| + |6|ln|6|-|a|(l-signasign6)ln|J|) , for |a| < |6| (190) 

(191) 



2 V ^ " ° ' '6 
iflln^ \d\ + pln^ |6| -aln|a| -61n|6| - if (|a|ln|a| + |6| In |6|) for \a\ - |6| < 1 



■asign6 1n|a| for |a| ^ |6| ~ 1 

For a; = 0, we see that the terms containing explicitly the complex factor i start always one power less in the loga- 
rithm compared to those without the factor i. The terms without the factor i start with the logarithmic contribution 
~ In Y^^-f-^, whereas the ones with the factor i start with terms ~ iJ^ In y^^^. Therefore, in the following section, 
where we apply the general formalism to the Kondo model, we will restrict ourselves only to those terms. This means 
that for ix) = we take into account all terms of order 0(1), 0( Jc) and 0{J^ In j^^^) without the factor i, and all 

terms of order 0{J^) and 0{J'^\Tij^p^^) with the factor i. Thus, we can neglect the contributions from L^^*^^*^ and 

2'(2c)c (|185p [according to Eqs. (|160p and (|16ip ]. and the first two terms on the r.h.s. of Eq. (|186p . Therefore, 

the first-order correction Z^^^ is not important for the final result. Furthermore, we can neglect all real terms in 
Eqs. p90p and (|19ip . This leads to the final equations 

Lf{E)t,=o = Lf'-EZ^^'^' , (192) 
Lf{E)^., « Lf^^-EZ^^^y-^lG^lr\^,,\G^l'-\G'^r^^^ (193) 



with 



Lf\E)t^=o « -iK f dxG^^^'Hx - i^) G^^^' Hx - i^) G'§' , (195) 

./n A^. A^. 



Ay = E,, - is , (196) 



and 



|6|ln|6| - |a|(l-signasign6)ln|f I for \a\ < \b\ 
-i I dxln{x-id)\n{x-ib) -i^ { |a|ln |a| + |6|ln|6| for jaj ~ |6J < 1 .(197) 

asign6 1n|a| for |a| <C |6| '--^ 1 








Close to resonance, where the effect of F can not be neglected, we have to consider the imaginary parts of the 
eigenvalues of Lg as well. If an eigenvalue h — iT occurs, we have to replace Ay — s- Eij — h + iT, and 

_ 2 jT^ . Ji 2 E h ~ 

I Ay I {Eij — h) — arctan , sign Ay — arctan^^^-= , ln|Ay| \n\Eij — h + iT] . (198) 



III. THE KONDO MODEL model in an external magnetic field ho > 0. In this case, 

we have 

A. Model and algebra in Liouville space 

Hs ^ hoS' , (199) 

Model. Now the RG formafism developed in the pre- _ 1 J ( J^^,)oS'V^g., for 77 = —77' = + /r,r,ri\ 

vious section is applied to the anisotropic spin-i Kondo ~ 2 1 — ( Jq/q)oS'V^,^ for 77 = —7/ = — ' ^ ' 
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where a denotes the reservoir index, i € {x, y, z}, and 

{JL')o = {JL')o = (^i-.')o . (201) 

S'* is the i component of the spin-^ operator of the quan- 
tum dot, (T* is a Pauh matrix, and {J^a')o [{Jaa')o] are 
the initial exchange coupUngs which correspond to pro- 
cesses without (with) spin flip. If we choose the reservoir 
states such that the exchange couplings are real, we find 
from hermiticity (|19p that the exchange coupling matri- 
ces are symmetric. 



(202) 



If one derives the Kondo model via a Schrieffer- Wolff 
transformation from an Anderson impurity model (see, 
e.g., Ref . |27[) . one finds 



= 1 



(203) 



which is used sometimes to simplify the calculations. Fur- 
thermore, although the general formalism and many of 
the following formulas are also valid for an arbitrary num- 
ber of reservoirs, we will finally consider the case of two 
reservoirs only, with a = L, i? = ±, and chemical poten- 
tials given by 



a- 



(204) 



where V is the applied voltage. 

Representation in Liouville space. In Liouville space, 
the initial Liouvillian and the initial vertex are given by 



r(0) 

'-'11' 



[Hs 
1 
2 



L 



for f] 
for T] 



(205) 



(206) 



where the spin superoperators LP — (L^^, L^^, L^^) are 
defined by their action on an arbitrary operator A in 
Hilbert space via 



L+A 



SA 



L-A 



-AS 



(207) 



We will now derive a closed set of basis superoperators 
to represent the renormalized Liouvillian and the renor- 
malized vertices (note that this set will be more complex 
than the one derived in Ref. [13 for the isotropic Kondo 
model) . Because the Hilbert space (spanned by the states 
I t) and I J,)) is two dimensional, the Liouville space of 
operators acting on it is four dimensional, and the su- 
peroperators defined in this section can be represented 
by 4 X 4 matrices. This means that we need at most 
4 X 4 = 16 basis superoperators. 

We define the four scalar superoperators 



(208) 



where 1 is the identity superoperator, and the vector 
superoperators 



= i (L+ -L- - 2iL+ X L- 



(209) 



Because the Hamiltonian is not isotropic, we have to split 
the vector superoperators (|209p into their components 
and consider these separately. We define 



L'^ = L"±iUy for i e {1,2,3,+,-} 



(210) 



and find that we get a convenient superoperator basis 
if we add the set L^^ L^, L^, L^} to the four 

scalar superoperators, where 



Li 



(211) 
(212) 



This set of 14 basis superoperators is sufficient to describe 
the anisotropic Kondo model with — — , where 
rotational invariance around the z axis holds. In the case 
of the fully anisotropic Kondo model with ^ J^, two 
additional superoperators 



(213) 
(214) 



have to be used, which we do not need in the following. 
We note the properties 



TrsL^ 



, for X = a, c, ft., 3z , (215) 
.± - , for x = 3,4,5 , (216) 
Trs , Tr5 L^' , Trs ^ , (217) 



TrsL? 



together with the transformation under the c transform 
(Ell) 

(218) 
(219) 
(220) 
(221) 

Using spin rotational invariance around the z axis and 
spin conservation, together with the properties (|215p - 
(PT7)) and Tt:sLs{E,uj) = [see Eq. we find that 

the Liouvillian and the current kernel can be represented 
as (in each step of the RG) 



{L^r = 




for X = 


a, 6, c, Iz, iz 










(Lir = 




for X = 


1,3 , 


{LIT = 




, for X 


= 4,5 . 



Ls{E,uj) = h{E,oj)L^' 
-iT''{E,uo)L'' - 



-iT''{E,uj)L'' 
iT^'-{E,uj)L^\ (222) 



^^{E,Lo) = iT\{E,oj)L'' + iT\\E,Lo)L^\ {222.) 

where as usual we always assume implicitly that the Trs 
is acting from the left when we consider the current kernel 
or the current vertex 



TrsS^ 



-'12 



Trs/:^ 



12 
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From Eqs. (flOTI) . (fT02| . ([STS]) . and (l2T9l) . we get the 

following transformation under complex conjugation: 

T^iE,uj)* ^ {- E, Lu) , ioT x = a, c,3z, {22A) 
T^E,uy = r^{-E,Lo) , iovx^b,lz, (225) 
h{E,Lu)* ^ h{-E,uj) . (226) 



The various terms in Eq. (j222p can be interpreted if 
one analyzes the spectral properties of the renormalized 
Liouvillian. hs{E,uj) has four eigenvalues, 

XoiE^Lj) = , 

Xi{E,Lj) = -iT''{E,ij) , (227) 
X±{E,uj) = ±h{E,uj) - iiV + r''){E,uj) . 

Ao corresponds to the stationary state, Ai describes 
the relaxation mode, and X± correspond to the 
two dephasing modes. The projectors Pi{E,u!) = 
\xi{E,uj){xi{E,uj)\ onto the four eigenspaces are given 
by 



Po(i?,a;) -L - j^a^j^^^^L , 
Pi{E,oj) = L° - + S^i^L^^ 



r°(i;,w) 



(228) 



1 



P±(i?,c.) = -(L^iL'^) , 



and the right and left eigenvectors follow from 

^^'■{E,uj) 
T'^iE.uj) 



{a(7'\xQ{E,Lu)) 


= Saa' 




{acT'\xi{E,w)) 


= Saa' 




'acj'\x±{E,Lo)) 


= <5.,- 




{xo{E,u:)\c7a') 


= Saa' 




{xi{E,uj)\a(T') 


= Saa' 




'x±{E,Lo)\aa') 


= 5a,- 


-a'Sa± 



(229) 



According to Eq. the eigenvector |xo(0, 0+)) corre- 
sponds to the stationary state. Therefore, we get 



1 



3z 



= -1 + 2MS^ , M = -— , (230) 

where M denotes the magnetization, which is related to 
T^^' (if no argument is written, we implicitly assume E — 
and cj = 0+). 

The stationary current follows from Eqs. (f45|) and (I223|) 

as 



which, by using Eq. ([230]) together with ^Trsi'"! 
TrsL^^S"^ = 1 and Tr5L''5^ = Tr^L^n = 0, gives 



To represent the vertices, we define the reservoir spin 
matrices 

= i(l±a^) , a± = ^ia,±iay) , (232) 

and the following operators in combined reservoir spin 
space and Liouville space of the dot 

~ , for X = a, 6, c, /i, Iz, 3z , (233) 

L± = i±crzp , for X = 1,3,4,5. (234) 

Using spin rotational invariance around the z axis, to- 
gether with the properties (fTOO| and (|2T5l) - (|2T7|) . the 
renormalized vertices can be represented as (for rji — 

-m = +) 

Gi2{E,UJ,UJi,Lil2)\ni=-r)2=+ — 

= E E G^\E,oj,oj,,uj2)L^ , (235) 

3z,3,4,5 

Ii2{E,uj,uJi,uj2)\-m=-v2=+ = 

= E E /^^'(i?,^,^i,^2)L? . (236) 

In these equations, we have used a compact matrix no- 
tation in the reservoir indices (ai,a2) and the reser- 
voir spin indices (cri,tT2). Whereas G^" = (G^^^^) 
jixs ^ (/^^aa) denote matrices in the reservoir indices, 
the quantities are matrices in the reservoir spin in- 
dices according to definitions (|233p and (|234p . Again, we 
note that Eq. (|236p holds only if the trace Trs is taken 
from the left. For rji = —772 = — , we use the antisymme- 
try properties ((57|) and and obtain from Eqs. (|235p 
and (ESS) 

Gl2{E,UJ,UJl,UJ2)\-qi=-V2 = - = 

= -E E G^'iE.u:,u;2,u:,f{L^f, (237) 

3s, 3, 4, 5 

/72(-B,t^,t^l,t^2)|,u = -')2=- = 

= -E E ^""^'(i?,^,^2,c.i)^(Lj)^ , (238) 

where (. . . )"^ denotes the transpose only with respect to 
the reservoir indices or the reservoir spin indices 

\^ Jaa' a'a ' V-* Jaa' a' a ' 



(LJ)L, = (LJ)<,v 



(239) 



(231) 



Using the properties p03p and (|104p together with 
Eqs. (|218p - (|22ip . we obtain the symmetry relations 

GZ'{E,LU,L0i,L02r = -G^f„(-i?,C^,C^2,C^l) 

for X = a, 5, c, Iz, 3z , 
G^i,{E,uj,uji,U2y = G^f„(-£;,a;,W2,^i), 

Gyj';'*(s, wi, W2)* = -gII,1;^''{-e,u,uj2,i^i) , 
G^(f;^(s, wi, W2)* = G^(^'~^(-S, w, W2, wi) , 

(240) 
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and the same for G ^ . 

Algebra. The results of a multipUcation of two of the 



operators Lj 



are summarized in the following table (all 
products not shown are zero): 



J 







fb 




fh 


flz 


f 3z 


L^ 






f 5 


LI 


LI 
u 




fb 


Li 
u 


Li 

U 




flz 


Li^ 
u 




f 1 

Lt 


u 


L% 
u 


u 


LI 
Ll 


Li 
Li 






Li 
Li 


Li 

Li 














L% 
TL% 






L% 
TL% 


Li' 


Li' 














Li 








±Li 





L^^ 













Li - Li 





±L% 











Li 
Li 


Li 





Li 


Li 



TLi 





±Li 









2L?;, 










Li 

Li 


Li 

Li 






Li 



TLi 









±Li 











L'^±L^ 


2L?^ - 2L^p 




(241) 



Note that, according to definitions (|233p and (|234p . the 
operators Lj_ are matrices in the reservoir spin indices, 
where each matrix element is a superoperator in the Li- 
ouville space of the dot. The same algebra holds without 
the reservoir spin indices if one replaces 



Li 
Li 



L^ 
Li 



for X = a, 6, c, h, Iz, 3z , 

for x= 1,3, 4, 5, (242) 



in Eq. (|24ip . with the only difference that also some 
products not shown are nonzero [these are the products 
— —LiLi defining the basis superoperators L±± 
according to Eqs. (|213p and (|214p . which are not needed 
for our case of rotational invariance around the z axis] . 
Furthermore, we note that the algebra (|24ip is also not 
changed if we replace 



Li 

Li 
± 



for X — a, 6, c, h, Iz, 3z , 

for X = 1,3,4,5, 
=F for all sign factors , (243) 



(Ll) 

-X 



which turns out to be very helpful to consider the two 
cases rj = ± for creation and annihilation operators, see 
the representations (|237p and (|238p . 

Finally, we note that the spin matrices (|232p fulfil the 
algebra (all products not shown are zero) 





^i 






^i 








^i 



(244) 



and the same holds if one replaces 



(245) 



B. Two-loop analysis 

1. RG above Ac 

Initial values. We start with the determination of 
the initial values of the Liouvillian and the vertices to- 
gether with their values after the first RG step, given by 
Eqs. (ESll-dZll). 

The initial values for the Liouvillian L^g^ and the ver- 
tices g[\] and g[\] follow from Eqs. pl5)) and pI5)) as 



r(0) _ 
Lg - 



[Hs,-]- = hoL'' , 



(246) 
(247) 



G« = Y.PG'i' = -ljoiL'' + L'^)a^ ,(248) 



where we have taken rj = —rj' = + and used a matrix 
notation Jq = [( Jg)ctQ'] for the exchange couplings in the 
reservoir indices. Inserting the various definitions of the 
basis superoperators, we find the following representa- 
tions for '>] = —'>]' = +: 



Ml) 
'-'11' 

^11' 



- JqsLJ — -Jo'iLt + LI) , 



(249) 



1 



Jo^Ll 



1 



Lf) + -J,iiLl + Ll), (250) 



where we sum implicitly over s = ± on the r.h.s. Ac- 
cording to Eq. ()4ip . the initial current vertex is given 



by I^IP = cJ^'Gii',. Using Eq. ^50\), we obtain for 



V 



-V' 



'ii' 



772 „ flz 



(251) 
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where 

(252) 

We have left out the components 3z and 3 in Eq. (|25ip 
because we assume imphcitly that the Trs acts from the 
left on the current vertex. 

To calculate the second-order contributions ((73)) - ([75l) 
from the first RG step, we use the following identities 
for any two superoperators Aiy = A^^LJ and Biy = 
B^'^L^, which can all be proven easily by summing over 
the two possibilities r]i = ±: 

An'B-ri = 2(Tr„AX^i?x'^') (Tr^L^Lj') , (253) 
^11' (£^11' — L^s^)Bi,i = 

^AZ'B^: (Tr.L?(i?o.' - hoL'^)Ly) , (254) 

with Eaa> = E + fXci — l^a'- Hcrc, Tr^ and Tr^ denote the 
trace with respect to the reservoir indices and the reser- 
voir spin indices, respectively. Applying these identities 
to Eqs. dZSll-jTll), inserting the results ([249l) - (|25T|l for 
the vertices, and applying the algebra of the basis super- 
operators, we find after some straightforward calculation 
the following result for the Liouvillian and the current 
kernel after the first RG step: 

2 

LUE) = -^^i2D + zE)(TTj^J^)L^ 

2 

-t^{2D + lE) (Tr(JoVJ - J^J^)) L= 

+h,{l ~ ^) (Tr JJJj) 

-i^ho (TtJ^J^) ^ (255) 

mE) - ^^S^-,{fIa - t^a') {{JL')' + 2iJL'f)oL' ■ 

(256) 

Thereby, the terms proportional to the band width can- 
cel with terms generated by the second RG step, see 
Eq. (|163p . Therefore, we can omit them together with 
the real terms of second order in J [which are not the 
leading order ones in this order, see the discussion before 
Eq. (|192p ]. Thus, we are left with the following inital 
condition to solve the RG equations ()158p and ()159ll : 

if)|A=Ao = -*|/io (XrJo^io^) (257) 

= cL.(A*a ~ Ma') {{JL'f + '2{jL'f),L' , 

(258) 

Z(2^)|a=a„ = Z^^''^U=A, = 0. (259) 

We now calculate the form of the vertex ([75]) after the 
first RG step, which leads to the initial condition for the 



vertices G^^?^^ and given by Eq. p45|) . We use 

the following identity, valid for rj = —rj' = +: 

AuB^v - (1 ^ 1') - 

^ ^x^^x'^' Ixlx _ ^x'^'^x^ ((i?)^ (L'^'fY -(260) 

Applying this property to Eq. (|145p . inserting Eqs. (|249p - 
(I25ip . and applying the algebra of the basis superoper- 
ators, we find after some straightforward calculation for 

T] = -rf = + 

^11' |A=Ao — ■^-'o '^O '^^s 

+ ^(A>Jo+Jo^Jo)Ll, (261) 

r7(2ai)| _ /"''^ P - P f'T^ -I- 

^11' |A=Ao — 8 "'0"'0 + 

+2JJ^ 4^ - 2 Jo^Jo^^)L^ . (262) 

Vertices in first order. The RG equations for the low- 
est order vertices G^\^ and are given by Eqs. (|115p 
and pT6)) with initial conditions ((249)) and ((25T|l . Using 
Eq. ()260p . one finds after some algebra that the initial 
form of the leading order vertices is preserved, i.e., for 
r] = —T]' = -I- 

g[\1 = ij^.L^-ij^(L^ + L^), (263) 

/7« = y-f^sLY + y-'^Ll, (264) 

with 

j:':/'- = ^i.'j:L^ , ^^,^-1(5,^-^), (265) 

provided that the renormalized exchange coupling matri- 
ces fulfil the poor man scaling RG equations 

(266) 

where I = InAg/A is the dimensionless fiow parameter 

and J^/^|a=Ao = Jq^'^ initial condition. 

The poor man scaling equations have the invariant 

(2cf = Tr(JV^-J^i^) = const . (267) 

They can be solved easily if one assumes the form (|203p 
for the initial couplings. In this case, we have the same 
form for the renormalized couplings, 

J'JJ- = 2^1-^ J'/^ , = 1 , (268) 

a 

and the two exchange couplings and fulfil the RG 
equations 

^J^ = 2(J^)2 , ^J^ = 2JV^ . (269) 
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These RG equations have two invariants 



Tk ^ A 



(270) 

where c = iV\iJ^7~~iJ^\ for \J^\ < \J-^\. Tk de- 
notes the Kondo temperature, which can also be written 
as 



K 



1 1 1,, 
A 6^57^ 55 



with 



5 = sign(J^) 



1- (V)' 



(271) 



(272) 



and, for |J^| < \J \ [where 6 = i sign(J^)|(5|], the 
imaginary part of the logarithm is defined such that 
— 27r < ln(z) < (such that Tk does not jump dur- 
ing the RG flow and is exponentially small). Thus, for 
I J^l < I J^l, we obtain 



where the flow parameter Z = In ^ = In ^ — In 7p- can 
be expressed in terms of the exchange couplings by using 
Eq. (j270ll : 



4c V + c 



ln> 



(281) 



Whereas the logarithms in this formula are pure 
anisotropy terms depending on the ratios of the exchange 
couplings, the prefactor ^ gives I ~ 0(1/ J), i.e., the sec- 
ond terms c^ Z on the r.h.s. of Eqs. p79p and (|280p are 
of the same order 0{J) as the first terms. 

Current kernel in first order. The current kernel 
T.^^ {E, Lu) can be determined similarly to L^^-* {E, lu), one 
just has to replace the first vertex on the r.h.s. of the RG 
equations (|130p and (|13ip by the current vertex (|25ip . 
Using the identities (|253p and (|254p together with the al- 
gebra of the basis superoperators, we obtain zero on the 
r.h.s. of the RG equations, i.e.. 



T, 



K 



Ke 



T (arctan | (5 1 + ^ {sigii( ) ~ 1 ) ) 



(273) 



sW(£;,c.) = 



(282) 



For the isotropic case (J^ = = J), we get the usual 
form Tk = Ae~2j. In terms of the invariants, the solu- 
tion of Eq. (|269p can be written as 



J" = c- 



1 



2c 



1- w 



(274) 



In the scaling limit Jq , Jq — > and Aq — > oo, such that 
4r = const and Tk — const (which is possible for > 
or I J^l < I J"*" I), we obtain the well-known isotropic form 

1 



r 



J' 



21n^ 



(275) 



with Tk given by Eq. (|27ip (where the anisotropy of the 
initial exchange couplings still enters). 



Liouvillian in first order. To determine L^g \E,uj) 

from Eq. (fT29l) . we use the solution ([263]) for g[]^, to- 
gether with the identities (|253p and (|254p to evaluate the 
r.h.s. of the RG equations and (fTM]) . We obtain 



Vertices in second order. The vertices in second or- 
der in J follow from the decomposition (I144p . where the 



first (imaginary) part iG^^?^'' (or il^y"''' for the current 
kernel) is given by Eq. (|145p . Since the renormalized ver- 
tices in first order have the same form as the initial ones, 
we obtain the same form as Eqs. (|26ip and (|262p with 



■7(2ai) 



J 



/=f(2ai) 



f7(2ai) 
^11' 



= -J^J^sL 



3z 



3 

s ; 



2i'^-Lj-L - 2J^J''^)L'l . 



(283) 



(284) 



In contrast, the (real) parts G^^^?^'' and I'^y"''^^ have 
to be calculated from the RG equation (I143P with zero 
initial condition. We use the same ansatz (I263P and (I264p 
as for the vertices in lowest order, i.e., for rj = —rj' = + 



dl 
d 



r(l) 



-{TrJ'J')hoL'' , (276) 
(Tr J-Lj-L)L" + 2c^L'' , (277) 



with the invariant c given by Eq. (|267p . Using Tr J^J^ = 
TiJ^J-^ + Ac'^ and J^J^ = fJ^ according to Eq. 
we find the solution (note that the inital conditions are 
zero in first order in J) 



4'\e^u;) 



(1) 



(E- 



(278) 



= iTr(J^ 

zw 



J^)hoL'' + 2c^ I hoL'' (?79) 



Tr(J^ - J^)L'' + 2c^ I L" 



(280) 



'-'"''^ - \k\sLl-\k^{Ll + Ll), (285) 



■7(2a2) 



2 "2 



(286) 



with - OiJ"^). However, instead of 

Eq. P55)) . we set 



0.0.' oo' oo' ' oo' 



(287) 



which can be viewed as a definition of B?'^l ^ . Inserting 
Eqs. (IMl), IMl), (EMU, (ESI, and (EHOl) into the RG 
equation (|143p , and using the following properties [analog 
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to ((2601) ]. valid for 77 = -77' = +: 

A12 B-^y - (1 ^ 1') = A^'^'b^"'" L^L^L^^', - 

-B^"'"A^''' (^{L^',f L^{L^','fY , (288) 
A23Bn'C-32 = 2BZ, (Tr^Ax'^'C^"^") • 

.(Tr.Lj'(LJ)..,ij';) , (289) 
I 



Together with Eq. (|287p . we have derived here the full 
two-loop equations for the vertex and the current vertex 
in the most general case, including anisotropics and all 
possibilities for the dependencies on the reservoir indices. 
If we define the total vertex by 

jz/± ^ jz/± ^ j^z/± ^ (294) 



we obtain after some lengthy but straighforward algebra 
the following RG equations: 



(290) 

(291) 
(292) 
(293) 



and take the sum of the RG equations (|266p . ()290|) . 
and (|29ip . we get, by neglecting terms of O(J^) and 
terms of 0( Jq J^, JqJ) (which vanish in the scaling limit 
Jo — > 0), the following result: 



Tr(J^ - J^) + 2c2 \{J"J^ + J^J') + \ (XrJ^J^) 

+\{k^j' + r'k^ + k''j^ + j^k'') , 
(TrJV") J^' - 2c2 I {r'^j^ + j^,r^) + w^j^ + j^iv^ , 



r 



-J'- = (l-Tr^^ 



{J' 



J' 



Tr(J- 



(295) 
(296) 



For ~ J — J and the case of two reservoirs with and 
Jll — JrRi these equations reduce to the two-loop equa- 
tions of Ref. [3^ [where the nonuniversal parameter a in 
Eq. (104) of this reference has to be chosen as a = 3]. 
In equilibrium (i.e., for a single reservoir), we recover the 
well-known two-loop RG equations for the anisotropic 
Kondo model derived in Ref.|36| [up to nonuniversal terms 
arising from adding terms proportional to the invariant 
Tr(» - J-^J-L)]. 

If the exchange couplings fulfil the relation ()203p . the 
RG equations can be solved easily. In this case, we have 

= 2V^^i^^/^ , (298) 

RYJ^ = ^1.'<^ , (297) R^^ = 2Vi^i?^/^ , (299) 
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with 





— 4(J — J 






+2(J^)3 4 


dl 


-4( - J 


1 r\ / T_L \ 3 




-2{r-f - 




dl^ 









(300) 
(301) 



+2{J^f + 2{J'K^ + J^K') , (302) 



(303) 



Special solutions of these equations are given by 

3 
2' 



+4/JJ(J-)^-(J-)^ln — 



(304) 



= I c^J^ + AlJ^rj^ - rj^ In — , (305) 

i?" = -{rf~2c^ir, R^ = -rj^ , (306) 

and the initial conditions K'^ = = R^ = R^ = can 
be fulfilled by adding linear combinations of the following 
solutions of the homogeneous part of the RG equations: 









(307) 


= r + 2i{j^f 




= j^ + 2irj^ 


, (308) 








(309) 



As a consequence, the final solutions are given by 
= ~ J^) + 2lc^r 



+lJ^{J^f-{J^f\n — 
+IJ^J'J^ - r-J^ In — 

R"- = - r-)-2c^ir- , 
R^ = j^{j§~r-) . 



(310) 



(311) 

(312) 
(313) 



Using Eq. ([^51]) and = (J^)^ _ (J^)'^, we see that 
the vertices in second order contain four different types 
of terms 

O(JoJ), O(J^), 0(J^ In J), 0(j2 In Jo) , (314) 

in addition to factors which contain ratios of J^/J^ or 
Jq/Jq- In the scaling hmit, the terms of order 0{JoJ) 
can be neglected. The terms of order 0{J^) are just a 
perturbative correction to the first-order vertices and can 
be neglected too. Therefore, we see that the terms R^/-^ 



are not important and the current vertex (j287p in second 
order can be written as 



K. 



cJ ,K 

OlOl' OLOL 



(315) 



The terms of order 0( In J) are logarithmic corrections, 
which are of order J In J smaller compared to the first- 
order vertices. Since J ^ 1 they are also perturbative 
corrections and can be omitted. They lead only to an 
overall change of the physical quantities without any in- 
teresting dependence on some physical energy scale (like 
e.g. the voltage, magnetic field, etc.). The most im- 
portant terms are those of order 0(J^ In Jo) since they 
diverge in the scaling limit. Therefore, they have to be in- 
corporated into the definition of the Kondo temperature. 
Using the solution (|274p with Tk T'k and expanding 

in J In TTT^ , we find 



(J^ 



4c 



(1 



= 2c- 



2c 



1 



J^ + 2(J^)^lng. 



J^ + 2JV^ In^ 



i.e., the terms of order 0(J^ In J^j'-) in Eqs. (|310p and 
(|31ip can be acccounted for by the redefinition 



Jq Tk 



(316) 



As a consequence, the main effect of the two-loop terms 
for the vertices is the replacement Tk T^, which we 
will always implicitly assume in the following. 

In other (more academic) cases, where the form (|203p 
of the exchange couplings is not fulfilled, the part R^^^J,^ 
can generally not be written in the form (|297p . i.e., 7 
must not necessarily be equal to a or a'. From the RG 
equations (|292p and (|293p . one can only prove that 



(317) 



However, in the special case of two reservoirs, these con- 
ditions lead again to the form (|297p . It is an open ques- 
tion of future research to analyze the two-loop equations 
in all cases and to find out whether there are interest- 
ing situations where the two-loop contributions can not 
be simply accounted for by a renormalization of a sin- 
gle parameter. In addition, it is not clear whether the 

corrections TC^^I'^ to the current vertex are generically 
unimportant. In this paper, we discuss only the physi- 
cally realizable situation, where the exchange couplings 
fulfil the property pI5)) . 

Liouvillian in second order. For Eq. (jl93p . we need the 
following second-order contribution to the Liouvillian: 



L 



(26) 



{E- 



(318) 
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where £5 •* and Z'^^''^ are determined by the RG equa- 
tions (|158p and (|159p with initial values given by 
Eqs. (P57)) and ((^55)) . Inserting the results and 
(pgg)) for the vertices into Eqs. ^55]) and ([15^)) . and us- 
ing the identities (|253p and (|254p . we find after some 
algebra the RG equations 

^4"'^ ^ -ITT ho (Tri^J^J") L^^ , (319) 

4^^"') = . (320) 
dl 

Using (J^)^ — Tv-^^' according to Eq. p66p together with 
the inital condition (|257p . we find the result 

Lf^ = -^^ho (TrJ^J^) L^^ , (321) 
Z^g''^ = , (322) 
where we have used that 

= Tr(JV" - J^J^) = Tr(JoVo" - Jo^ Jo^) 

is an invariant according to Eq. ()267|) . 

Current kernel in second order. Analog to L^g'^\E,uj), 
we evaluate the current kernel 

{E, u;) _ ^ (323) 

from the RG equations (|158p and (|159p by using the re- 
sults (|264p and (|284p for the current vertices. This gives 
the RG equations 

= _,|cL.(A^a-A^a'){^^.'(^^^^)«'« + 

+ J^„,{r-J^ + J^r')a'a}L' ,(324) 

and ^4^''' = 0. Using J^J^ ^ iJ^ ^nd JV-^ -|- 
= 2^ J"*- according to Eq. (|266p . we find together 
with the initial conditions (|258p and (I259P the result 

) = JcL,(Ma - M.') [(J^.')' + 2(Ji-„0'] , 

= 0. 

(325) 

2. RG below 

We now start to evaluate the final formulas p92p - (|195p 
for the effective Liouvillian and the current kernel, which 
are represented in terms of the basis superoperators ac- 
cording to the forms (|222p and (|223p . All final quantities 
are evaluated at w = and are meant at A = which is 
not indicated explicitly, i.e., we write from now on 

LsiE) = Ls{E,uj ^0)a^o = 

= h{E)L^' ~iV{E)L'' 

~ir''{E)L'' -ir^'{E)L^\ (326) 
^^{E) EE S^(£;,co = 0)a=o- 

= iT''jE)L'' + iTl^''{E)L^^ . (327) 



Since only the first-order vertices evaluated at A = Ac 
occur in the final formulas, we will furthermore use the 
convention 

Jaa' = (•^qq')|a=A, , Jaa' = (•^qq')|a=Ac ■ (328) 

According to Eq. (jllip . the scale Ac is defined as the 
maximum of all physical energy scales, i.e., in our case of 
the Kondo problem 

Ac = max{i;, U, h} , (329) 

where h is defined by the renormalized magnetic field or, 
more precisely, as the real part of the eigenvalues of the 
operator Ls- Therefore, we have to keep in mind that the 
renormalized exchange couplings p28p depend implicitly 
on the variables E, V, and h via the scale Ac. 

According to Eq. ([55]) . any function H{Ai2) of the 
quantities A12 = £'12 — Ls is defined by 

H(Ai2) - ^H(£;i2-z,)P,(2^) , (330) 

i 

with i = 0,1, ± and Zi — Xi{zi). Xi{z) = \{E,uj) and 
P,{z) EE P.i{E,uj) are given by Eqs. ([#7)1 and with 
z = E + iLu. As a consequence, we obtain the following 
self-consistent equations for z^: 

^0 = 0, (331) 
zi = -iT^zi) , (332) 
z± = ±h{z±)-i{r'' + T^)iz±) . (333) 

The pole zq = corresponds to the stationary state, 
which does not occur in the resolvents between the renor- 
malized vertices due to our special construction, where 
the eigenvalue zero is perturbatively integrated out dur- 
ing the first discrete RG step. The poles zi and z± corre- 
spond to the spin relaxation and dephasing modes. Using 
the symmetry properties (|224p and (|226p . we obtain un- 
der complex conjugation 

z* = ir"(-zi*) , (334) 
zl = ±ft(~z;)+z(r'^ + n(-zi) , (335) 

which gives 

zi = —z\ = — iFi , (336) 
z± = -z; ^±h- if2 . (337) 

As a consequence, Eq. (I330p becomes 

H(Ai2) = 7^(£;i2)Po(0)+7^(i;i2+«fi)Pi(-ifi) 
+n{E^i~ sh^iV2)P,{sh~%V2). (338) 

where we sum over s = ±. 

Inserting Eq. (j338p for the various functions occurring 
in Eq. (|193p - (|195p . we see that the projectors stand al- 
ways left to G^l^" or cf"'^'^''. Using the results (gMl) and 
(I283P for these vertices together with the algebra of the 
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basis superoperators, we conclude that the parts of the 
projectors (|228p containing either or L^^ do not con- 
tribute, and we can replace ^ 1. Thus, by inserting 
Eq. in Eq. we obtain the useful identity 



components of the effective Liouvillian in second order 



W(Ai2) ^ 



H{Ei2 

- sh 



(339) 



which will be used frequently during the following anal- 
ysis. As a consequence, we find that the spin relaxation 
rate Fi will cut off all logarithms where the magnetic field 
does not occur, and the spin dephasing rates T2 cuts off 
all logarithms where the magnetic field occurs. This is 
expected since the spin dephasing rate corresponds to 
spin flip processes. In addition, we also conclude from 
Eq. I|339p that the pole of the projector Pi{E,lj) at 
r°-{E,uj) = does not contribute. 

Liouvillian and current kernel in first order. The Li- 
ouvillian and the current kernel in first order in J have 
already been determined in Eqs. (P75)l - P5D)) and ([^5^ . 
which gives 



h'-'HE) 
V'-^^E) 

r<^\E) 



^TT{r - J^)ho + 2c''l,ho, (340) 



-2iEc^L , 







(341) 
(342) 
(343) 



with /c = In^, and 



r^(i)(^) - tI^'^-^\e) = . 



(344) 



Liouvillian and current kernel in second order. The 
first two terms on the r.h.s. of Eq. (|193p have already 
been evaluated in Eqs. ([5^ . ([5^ . and ((5^ . For the 
evaluation of the last two terms on the r.h.s., we use 
Eq. (|339p and the identity [Ti-iE) is any function] 



n{En'+z)An' Byi = 

= 2n{E^^, + z) Tr.Lj L^Lj: , (345) 



together with the results (|263p and (|264p for the first- 
order vertices. After some algebra, we obtain for the 



Re/i(2)(£;) = 

Ivah'^'^XE) = 
ReF'^(^)(£;) = 
ImF'^(2)(^) ^ 



+ {E -~E) , (346) 
-h\{j:^,f-{E^-E),{M7) 

^\E^a' - h\{jL' f + iE^ -E) , (348) 



'-\E, 



-{E~.^E) , 
Re[F"(2)(i;) +r"(2)(^)] ^ 

^l{\Eo.,.. ^h\{J^^,f + \E^,..\{jL'f] 

+ {E^-E) , 
Im[F''(2)(£;) +F'=(2)(i;)] = 

= -J {{Eac.' - h)C2{Ec.a' - h)iJL' f + 
+Eaa' ^l{Eaa')iJaa'y} — (E ^ —E) , 

r^^(^){E) = ^hoTTj^J^ , 



(349) 



(350) 



(351) 
(352) 



and the following result for the components of the current 
kernel in second order: 

rf \E) = - J(^„ - ^^^,) [JZ,JL, + 2Jli,jt^) , 

(353) 

ReVl^^^^KE) = ^\E^^. - h\Jlt.JL' + {E^ -E) , 



(354) 

h\C2{E^o.'-h)JZt,JL' 
^{E^~E) , (355) 

where we have used the shorthand notation 

An 



A(x) 



In 



(356) 



for the logarithmic terms. Note that the results respect 
the symmetry properties (|224p - ()226p . i.e., the real parts 
of all quantities are symmetric in the Laplace variable E, 
whereas the imaginary parts are antisymmetric in E. 

As we can see, for finite Laplace variable E, logarithmic 
corrections occur already in second order in J. This is 
the case for the real part of h^'^\E) and the imaginary 

parts of V<^\E), V'^'^^E), and tI,''^'^\e). Therefore, 
in the stationary state (i.e., for = 0) no logarithmic 
contributions occur in 0{J^) for the conductance and the 
magnetization, whereas the time evolution is influenced 
by logarithmic contributions already in this order, see 
Sec. lIV Al for more details. 

Liouvillian and current kernel in third order. In third 
order, there are two contributions to the effective Liou- 
villian and the current kernel, denoted by the superscript 
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(3a) and (36) and given by the two expressions, Eqs. (|194[) 
and (|195p . The components (3a) from Eq. (|194p can 
be evaluated simihar to the components in second order 
from Eq. p93p shown above, just by using in addition 
the results (|283p and (|284[) for the imaginary part of the 
vertices in second order. We obtain after some algebra 
for the components of the effective Liouvillian 

/j(3a)(£;) ^ T''^^^'> (E) = T^^'^He) ^ , (357) 

p3.(3a)(^) = -^{E^c.' -h)MEac.' ~h)■ 
■JL'iJ"J^ + J^J'U'c + {E^ -E) , (358) 
and for the components for the current kernel 

Yf-)^E) = -'^[e^o.'Ci{E^o.')JZ,{J^J^)c.'o.+ 

+ {E-.-E) . (359) 

As expected, the logarithmic terms for the real parts of 
F'^^ and start in order . This means that the loga- 
rithmic terms for the stationary magnetization (|230|) and 
the stationary current (|23ip start in third order in the ex- 
change couplings. Note that we have not calculated the 
terms in order 0{J^) without a logarithmic contribution. 

The remaining logarithmic contributions in third order 
are contained in the component (36) , given by (|195p and 



((T97)) . With A = G(i)^ and B = Jt(i)^ we have to eval- 
uate expressions of the form [note that {A^'^Y' — ^ ^'^d 
(B^")'^ ^ -B according to Eqs. ([263]) . ([264ll . and (|202| ] 

(^/B)l2i?(Ai2)A235(Ai3)^3l = 

•Tr,{LXi?(A,^„jLj;5(A„,„3)L^;;- 
-(Lj)^i?(A„,„J(Ljy^(A„3„J(Lj:)^} , (360) 

where A^.q. — Ea^a - — Ls, and the two functions R{E) 
and S{E) are either both symmetric or antisymmetric: 
R{-E)S{-E) = R{E)S{E). This identity can be de- 
rived by summing over the two possibilities rji = ± and 
using the representations ([235| . (|236| . (|237| and ([238| . 
The functions {R/ S){Aa^aj) on the r.h.s. can be ex- 
pressed in terms of the basis superoperators by using 
Eq. ((3391) . 

(i?/5)(A„^„J ^ {R/S)^,^YE + tr,){l-L-) 

+ {R/S)a.,a., {E-Sh + zf2)i(L^ + sL'') , 

(361) 

where Raioj {E) — R{E-\-^ai —fJ-aj)- Inserting this result 
in Eq. p60p . we find after some straightforward algebra 



G^f i?(Ai2)G^f 5(Ai3)G|^f = -^J^,o.2JtaXso..Rc.,c.2{s'E - sh)S^,^,{s'E)[L^ -f ss' L'^) 

-\j^,o.X.o.,Jl,c.,Ro.,o.As'E)So.,^,{s'E - sh){L- + ss'L'^) 
-^^io.^a.c.3'^i3a,^c.,..(s'i?-s/i)5„,„3(s'i?-s/i)(L'^-L^) , (362) 



=.(l)c _ 



1 



and for i5 = 



/7f)^i?(Ai2)G^;^^5(Ai3)G^f l^^=o = ^ic2;,„,J:,„,J„^,„3J„^3„,i?o,„,(0)5„,„3(s/i)Li 



^ 2 ''0102 "^0102 '^0203 "^Q3Q1 ^aia2 (s^)'5'qiQ;3 {sh)L 



(363) 



Using this result in Eq. (|195p . we have to distinguish 
many different cases in order to evaluate the logarithmic 
contributions according to Eq. (|197p . To simplify the dis- 
cussion, we will assume in the following only two reser- 
voirs with chemical potentials = Q:y. Furthermore, 
we take w.l.o.g. V,h > 0. As we will see in Sec. IIV Al we 
need only the two cases E — and E — h in Eq. (|362p 
for the calculation of the renormalized g factor and the 
spin relaxation and dephasing rates up to the first 
logarithmic correction. For the magnetization and the 



current, we consider the stationary case E ^ 0. 

We start with E = 0. In this case, there are four 
possibilities for the parameters \a\ and |6| in (|197p : 



- V h \V±h\ 
1""'"^" = AT' Al'^V— 



(364) 



According to Eq. (|197p . logarithmic contributions can 
only occur if at least one of these parameters is small 
compared to one. Therefore, we can distinguish three 
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cases 



V <$:h 



Ac 

\V-h\ 



Ac 

IT < 1 

Ac 



< 1 



(365) 



and all the other parameters not indicated are of order 
one. Discussing these three cases separately and collect- 
ing only the logarithmic terms according to Eq. (|197p . we 
find after some lengthy calculation 



± \2 



+ -\V ~ h\C2iV ~ h)J^iJ^^) 



pa(36)^Q^ 
TT 



2 



+-{y-h + \v- h\)C2{v - h)j^^^Jt^Jt . 



(366) 



(367) 



-\\v-h\c,{v-h),rJjZjl 



+-{y-h)C2{V-h){JZJt^ 



(368) 



with 



jz/± _ jz/± jz/± _ jz/± _ jZ/± 



'lr 



'LR 
'RL 



'RL 



(369) 

(370) 



All other components do not contain any logarithmic con- 
tribution in third order in J. 

Next we consider E = h. In this case, there are six 
possibilities for the parameters \a\ and |6| in (|197p . 



- V h \V±h\ \V±2h\ 



Ac ' Ac 



Ac 



Ac 



(371) 



In addition to the three cases shown in Eq. ()365|) . we 
have to consider the additional case 



Liouvillian aX E = h 

n ~ 



r 



±\2 



(373) 



-^6{h -V){V~ h)C,{V - h)J^^iJZr, 
-\b{V ~h)iV~ h)C,{V ~ h)Jl{Jtaf ■ 



IV. RESULTS 



(374) 



(375) 



In this section, we summarize the zero temperature 
results for the anisotropic nonequilibrium Kondo model 
in a finite magnetic field, introduced at the beginning 
of Sec. IIII Al We will calculate quantities characteriz- 
ing the exponential time decay of the magnetization, i.e., 
the spin relaxation and dephasing rates T1/2 together 
with the renormalized magnetic field h. Furthermore, we 
present results for the stationary current (/''')''* and the 
stationary magnetization M. All quantities are calcu- 
lated one order beyond leading order, i.e., up to the first 
term leading to logarithmic enhancements (suppressions) 
at resonance. 

For simplicity, we evaluate the results for two reservoirs 
with voltages 



Mo 



V 



L/R= ± 



(376) 



The diagonal and nondiagonal exchange couplings are 
denoted by 



jz/± ^ jz/J 



J, 



z/± 



J 



LR 



J 



z/1^ 
RL 



(377) 



The various couplings are renormalized couplings eval- 
uated at (note that the Laplace variable E is fixed to 
either = or = /i in the following) the cutoff scale 



Ac = max{y, K] 



(378) 



F~ 2/1 



\V-2h\ 
Ac 



< 1 



(372) 



Discussing these four cases, we find after a lengthy cal- 
culation the following logarithmic terms for the effective 



Unrcnormalized (bare) couplings are denoted by the in- 
dex "0". Furthermore, we consider w.l.o.g. the case 
V,h > 0. Although our results of the previous section 
include all cases for the ratios between diagonal and non- 
diagonal couplings, we will treat in this section the real- 
istic case where the exchange couplings fulfil the relation 
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(I203)) . i.e., 



XL+XR = 1 . (379) 



In this case, the two couphngs and can be calcu- 
lated analytically from Eq. (|274p with 

(380) 

according to Eqs. (I270p and pi6p . Furthermore, we note 
that some of the third-order terms can be taken together 
because 



TZ ( J±\2 _ 7^ 72 7^ 



(381) 



Finally, for the logarithmic terms we use the shorthand 
notation 



Ci{x) — In 



A, 



L2{x) = In 



A. 



x^ + Tl 



(382) 

and close to resonance, we use the broadened sign func- 
tion 



• / N 2 X 
sign,- X — — arctan — 



(383) 



where f i (f 2) has to be chosen when the magnetic field 
does not (does) occur in x (this rule refers to the case 
E = {), ior E — h also other cases can occur, see below). 
Furthermore, we use the broadening of the sign function 
also for expressions involving the absolute value or the 
theta function via 



\x\^ = xsigHiix) , 
^iix) = -[1 + sign^(a;) 



(384) 
(385) 



A. g factor, spin relaxation and dephasing rates 

In this section, we determine the renormalized mag- 
netic field h or the renormalized g factor 



= 2- 



dh 



(386) 



together with the renormalized spin relaxation and de- 
phasing rates Ti/2. These quantities enter all formulas, 
especially they determine the resonance position V — h 
and how the logarithmic divergencies are cut off by the 
rates. We note that previous works have only calculated 
ri/2 up to 0{J^) in bare perturbation theoryj-^ and the 
renormalized magnetic field has so far been calculated 
only up to 0{J), see, e.g., Ref. [H. In contrast, here 
we perform a renormalized perturbation theory for the 
rates and calculate in addition the logarithmic terms of 
0{J^ In) for Ti/2 and those of 0{J^ In) for h (correspond- 
ing to contributions in three-loop by using the conven- 
tional classification). Furthermore, we emphasize that 



we do not need to combine self-energy terms with ver- 
tex corrections as in slave-particle formalism^ since we 
have directly set up a kinetic equation from the very be- 
ginning where the renormalized rates and the renormal- 
ized magnetic field can be directly read off by studying 
the poles of the reduced density matrix of the quantum 
dot in Laplace space, see the discussion at the end of 
Sec. Ill Al We now calculate h and F1/2 perturbatively in 
the renormalized exchange couplings by using our results 

(iMUD-dM^. dsinii-iissii), dssi), (isn3, and mm-mm 

for h{z) and r°/'^(z) (note that we take the analytic con- 
tinuation E E + iLU oi these equations). Thereby, we 
consider all terms of 0{J^, In) for ri/2, and all terms 
of 0(JO), 0(J), 0(J2 In) for h. Terms of 0{J^) [0{J'^)] 
without any logarithmic contribution arc consistently ne- 
glected for ri/2 (h). 

To get the perturbative solution, we insert first the low- 
est order results into Eqs. and ([555)1 . 
This gives 



Zi — —I 



1- J5 + (J^)o 



(387) 



1 + UJ^-{J^)o) + 2c\ 



ho 



+ (/»^''H2+)-«(r"^''' +r^^''^)(2+)) . (388) 



Using the forms (|336p and (|337p for zi and z+, we see 
that ri/2 starts at 0{J^). Therefore, neglecting terms 
of O(J^), we can replace zi —> and z+ —> h in the 
arguments of F"/'^^*^^ and /i^*^^ for k = 2,3. Furthermore, 
we can neglect all contributions of 0{J^) for the first 
term on the r.h.s. of Eq. (|387p . and all terms of 0{J^) 
for the first term on the r.h.s. of Eq. (|388p (since they 
do not contain any logarithmic contribution). Thus, we 
finally obtain 

fi = F''(2)(o)+r°(3)(o), (389) 
T2 = J2 {-im/i(''HM + Rc(r''('^) +F'=('^))(/i)} , 



fe=2,3 



(390) 



where 



(392) 



is the renormalized magnetic field up to first order in J. 

Inserting the resuffs (p46)) - ((35T|) . ((366)) . ((367)) . and 
(l575 )) - ((575)) into Eqs. (P5g)) - (|5gT|) . and specializing to the 
case (|376p of two reservoirs and the special form (|379p of 
the exchange couplings, we obtain 
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fi = lh{jif + ^-{\V -h\, + V + h){j^^f 

WhC2{h)r^ {j^f + l\v- h\,C2iv - h)j^ - - h)C2(y - h)J^^J^^J^ , (393) 
f2 = ^v{j:i^f + ^h{j^)' + ^[\v-h\, + v + h) (Ji,)' 

(394) 

h = h- hc^{h) {J^f + h)C,{V - h) {Ji^f . (395) 



These formulas show precisely which rate cuts off the 
various logarithmic terms. Since the determination of the 
renormalized rates and the renormalized magnetic field 
involves also the value of the Liouvillian Ls{E) a,t E — h, 
we can not just use the simple rule that Fi (r2) cuts off 
the logarithms where the magnetic field does not (does) 
occur (as it is the case for the stationary case E — 0). In 
contrast, we find that the logarithmic terms of h and r2 
are cut off by Fi, and those of Fi by F2. 

h is the well-known renormalized magnetic field up to 
first order in J which determines the precise resonance 
positions. Note that the sign of the linear terms in J in 
Eq. p92)) is different from the one in Eq. (|279[) . show- 
ing that it is important to consider the linear term in 
frequency of the renormalized Liouvillian to determine 
the correct renormalization of the magnetic field (sim- 
iliar to field-theoretical rescaling techniques from the Z 
factor-^"*). Note that the renormalized magnetic field is 
linear in the renormalized coupling, but starts only in 
second order with respect to the bare coupling. Further- 
more, we note that the determination of the renormalized 
magnetic field up to first order in J is usually associated 
with a two-loop calculation, see, e.g., Ref.lssj i.e., the first 
diagram of Fig. 2] is called a two-loop diagram. We are 
not adopting this notation in the present paper because 
this diagram arises from closing the one-loop diagram for 

I 

dho 
dh 

dho 
~ „ dh 




where we have used Eq. (|38ip to take some terms to- 
gether. 

For both rates Fi/2 5 we get a jump in -gj;^ a,t h = V 



I 

the renormalization of the vertex (the third diagram in 
Fig. 2]). Therefore we call this diagram also one- loop and 
the second diagram of Fig. 0] a two-loop diagram (which 
usually would be classified as three- loop). We note that 
the classification in one-loop and two-loop contributions 
is not unambiguous, e.g., if one uses a not-normaled or- 
dered version of the RG formalism, the renormalization 
of the magnetic field in lowest order in J would arise 
from closing the renormalized vertex with itself, i.e., in 
this case one would classify it as one-loop. Therefore, 
we prefer in this paper a classification according to the 
orders in J and not with respect to the loop topology of 
the diagrams. 

Some insight into the logarithmic terms in Fx/2 ^^^d 
h can be gained by differentiating these quantities with 
respect to the bare magnetic field Hq [for the rates F1/2 
and the second-order terms in h, it is actually equiva- 
lent to differentiate with respect to h because multiply- 
ing with ^ 0{J) would just bring the expression one 
order higher in J, leading to terms which we have ne- 
glected anyhow]. Furthermore, we differentiate only the 
prefactor of the logarithmic terms and disregard the de- 
pendence of the renormalized exchange couplings on 
(which are small corrections since the logarithm varies 
only slowly). This gives 



(396) 
(397) 
(398) 

I 

in the leading order. In the next to leading order, there 
is a logarithmic enhancement for — > for both rates 
and a logarithmic enhancement for h > V in while 



= f {J^Y + -^Hh - V) {J^^f + 7:C2ih)J^ {J^f + ne2{h - V)C2{V - h)J^ (J„^d)' , 
= 2-{Ji-{rM)-Cm{jtf-L,{V-h){jt^f , 
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0.005 



I 1 1 1 1 

0.5 1 1.5 2 

ho/V 

FIG. 6: The relaxation and dephasing rates Fi and de- 
rived with respect to the magnetic field ho, for the isotropic 
Kondo model with V = IQ'^D and Tk = IQ'^D. ^ ex- 
hibits a logarithmic enhancement for h > V whereas is 
suppressed for h < V. 



1.95 




1.91 I 1 ' 1 

0.5 1 1.5 

ho/V 

FIG. 7: g factor g — 2dh/dho, derived with respect to the 
magnetic field ho, for the isotropic Kondo model with V — 
IQ-'^D and Tk = 10"*D. 



we find a logarithmic suppression in for h < V , see 
Fig. El 

The renormalized g factor shows a logarithmic sup- 
pression for h ^ Q and another logarithmic suppression 
for h V which is symmetric, in contrast to the ef- 
fects we observed for the rates, see Fig. [71 As one can 
see, the second-order logarithmic terms lead to signifi- 
cant changes of the first-order shift of the g factor. Es- 
pecially at resonance V — h, where one usually reads 
off the renormalized magnetic field from the conductance 
(see below), the change in the g factor by the logarith- 
mic contributions is very important. However, fixing the 



voltage d.iV = h, we obtain 

(399) 

with 

{^i)\v=H^ly[{Jif + %Jt^'] ■ (400) 

As a consequence, the g factor aX V = h shows only a 
weak dependence on the bare magnetic field /iq via the 
renormalization of the exchange couplings. The second- 
order terms just lead to an additional overall decrease 
^ ,P In J. To see the suppression of the g factor more 
clearly at resonance, it is necessary to measure it di- 
rectly for various values of ho/V. In principle this can 
be achieved by electron spin resonance (ESR), but, si- 
miliar to previous measurements of the splitting of the 
Kondo resonance in the spectral density^^ we propose 
here a simpler setup with a weakly coupled third lead, 
see Fig. |8l We use the form (|203l) and, for simplicity, 
choose a symmetric coupling to the left and right leads, 

XL ^ xr ^ X , < 1 , (401) 

'^P,nd — "^PL ^ ^PR ■ l^U^J 




Jlr 



FIG. 8: Three-terminal setup including a third weakly cou- 
pled probe lead, with fip = eVp, in order to measure the 
renormalized g factor as function of the applied voltage eV = 

By changing the voltage Vp of the probe lead and mea- 
suring the differential conductance Gp — one can 
probe the renormalized g factor as function of the voltage 
V = — fJ-R between the strongly coupled left and right 
leads. The reason is that the contribution of inelastic co- 
tunneling to the probe current changes when the probe 
voltage crosses the resonance points Vp = i /i, — y 
The probe current up to second order in Jp^nd can be cal- 
culated from Eqs. (p3T|) . ([353ll . and ([SM]) with the result 

(/^)«* = rp('^(0) + 2Mrp('^(0) , (403) 
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where 

^6(2 



(0) 



Vp 



r 



lz{2). 



V 



(404) 



V 



(0) = _ _ h\2 + \Vp + --h\ 



\Vp-'^+h\2-\Vp + ^+h\2 



'P.nd 



(405) 



Here, M is the magnetization, which in lowest order is in- 
dependent of the probe lead and is calculated in the next 
section. Taking the derivative with respect to Vp, we 

find that '^^^y^'^^ is independent of Vp , whereas ^'""^y 
gives rise to a steplike structure shown in Figs. [51an(f [TUl 
for the two regimes V > 2h and V < 2h, respectively. 
Sufficiently far away from resonance, we obtain analyti- 
cally from Eq. ([405]) for V > 2h 



dVp 



h>Vp> ^-h 



p.nd) 




and for V < 2h 



2, for 



v_ 

2 



h>Vp> 



1, ioY ^■ + h>Vp> -^-^h 
or ^-h>Vp> -^-h 



0, otherwise. 



(407) 

From the width of the steplike features one can read off 
the renormalized magnetic field h at given voltage V and 
magnetic field h^. Thus, by varying V or /iq, one can 
measure the renormalized g factor as function of ho/V, 
thereby revealing the suppression at resonance V = h as 
shown in Fig. [71 We note that the experimental situa- 
tion is usually in the strong coupling regime J ~ 0(1), 
where the qualitative features are expected to be more 
pronounced. 



B. Magnetization and susceptibility 

The stationary magnetization (which is nonzero only if 
a magnetic field /iq 7^ is applied) can be calculated from 
the stationary density matrix (|230p . This density matrix 
is the solution of the kinetic equation and is a diagonal 
matrix. The occupation probabilities for the states | t) 
and I) are 



2r"(0) ' 2r^{0) 

and the magnetization is given by 

1 , , r3^(0) 

M^-(PT-n) = -^ . 



r'^(o) + r3^(o) 



(408) 
(409) 



0.022 



0.02 



0.018 



V >2h 









*- 2h -*• 


<*-V-2h-P- 


4- 2/i 









0.016 

-5 -4 -3 -2 -1 1 2 3 4 5 

FIG. 9: Probe-lead conductance in a three-terminal setup for 
ho = 10"*D, V = 5 X IQ-^D, and Tk = 10"*D. 

0.035 



V <2h 

0.03 
0.025 
0.02 
0.015 
0.01 
0.005 

-2-1012 

fip/ho 

FIG. 10: Probe-lead conductance in a three-terminal setup 
for ho = IQ-'^D, y = 8 X IQ-^D, and Tk = 10"*D. 



V 



-2h - V- 



V 



r°(0) = fi is given by Eq. ((3931) . r3^(0) can be calcu- 
lated from Eqs. ([55^ and ([5551) : 



(0) - lh[{j^f + 2{jt,f} 
WhC2{h) { {Jif + J^K^JL) (410) 

-|(y - h)C2{v - h) { {Jt^f + JtJl^JL) , 



where we have replaced hg h in the first term, which 
gives only rise to negligible terms of 0(7^,7"* In). Be- 
cause we have considered all terms of 0(J^, J"^ In) for the 
rates, we can calculate the leading order 0(J") (which 
is independent of the couplings) and logarithmically en- 
hanced terms in O(Jln) of the magnetization. 
Using Eq. (I38ip and the abbreviation 

X = HJ^f + 2h{Jt^f + 2hC2{h)Jl{J^f 

-2{V-h)C2{V-h)J^(J,i,)\ (411) 
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we can write the magnetization in the compact form 



M 



(412) 



X + 2hC2{h)Jl{Ji 



± \2 



X^{\V -h\^ + V -h){l + C2{V 



We note that, in the isotropic case, and when expanded 
systematically in J, the terms of 0(J°, Jin) are in agree- 
ment with the perturbative results obtained in Ref. [13 
and the poor man scaling results of Ref. [13 if the differ- 
ent definition of the magnetization is taken into account 
(the magnetization in those references is M = 1 in the 
ground state, whereas it is M — —1/2 here). However, 
we have shown here how the spin dephasing rate cuts off 
all the logarithmic divergencies (the treatment in Ref. [g^ 
was only perturbative in the bare coupling and did not 
treat the case of finite magnetic field). 

If the voltage is smaller than the renormalized mag- 
netic field /i, the logarithmic term ^ £2(^1) in the nu- 
merator of Eq. (|412p must be consistently neglected and 
the magnetization is equal to —1/2. This is because the 
energy difference provided by the voltage is insufficient 
to flip the spin 1/2 in the dot out of its ground state | I). 

In the limit V ^ h, the logarithmic term ^ C2{V — h) 
in the denominator of Eq. (|412p must be neglected and 
the magnetization is given by 



M 



which, for J^^^ ~ Ja ^ ^ can be expanded as 



1 x + 2hC2{h)ri,{Jt^f 

2 X + 2V{Jt^Y 



(413) 




+C2{h).P^\l+[^ 




(414) 



As we see, there is an interesting logarithmic enhance- 
ment at /i = 0, which is a pure nonequilibrium effect [if 
we set = in Eq. (1413]) . we get M = -\]. 

In the limit V > h, V — h <^ h, (i.e., for voltages which 
are slightly larger than the renormalized magnetic field), 
we can neglect all logarithmic terms ~ C2{h) and find 
after expanding in J 



which we calculate by differentiating the magnetiza- 
tion (|409p with respect to the bare magnetic field ho 
[or the renormalized field h by multiplying with = 
1 + 0(J,J^ In)]. In the special case of symmetric cou- 



plings {J^ 



•^nd ^'-'^ X ~ -2, -L), the magnetization 



foT V > h can be expanded as 



M = -- 



V+h V+h 
-C2{h)J'' 



^ lc2{V~h)r 



[v + hY 

and the susceptibility becomes 

V 2V 



(417) 



X 



2V{V -h) 



C2iV~h).F 
C2{h)J' . (418) 



To visualize the logarithmic effects better, we multiply x 
with ~{V + h)^/V and get 

(V + hf^ ^ ^ ^ 2rC2{V -h) + 2^^—t,rC2ih). 



V 



(419) 

Plotting Eq. (I419p as function of the magnetic field ho 
emphasizes the logarithmic enhancements at /i <C and 
h^V, see Fig. [TTl In this plot, we multiply the suscep- 



X 




M 



V-h 



l + iJl+J^)C2{V~h) 



(415) 



This result shows the logarithmic enhancement of the 
magnetization slightly above the resonance. 

From the magnetization, we can also calculate the sus- 
ceptibility 



X 



dM 
dho 



dM 
'dh 



(416) 



FIG. 11: The susceptibility x, multiplied with , 
as function of the magnetic field ho for the isotropic Kondo 
model with V = lO""*!) and Tk = 10~**D. 



{V+haY 



because 



tibility x not with —- — but with y 
this change causes only subleading corrections and be- 
cause the unrenormalized magnetic field ho is easier to 
access experimentally. 

In summary, from the various expansions Eqs. (|414p . 
(|415p . and (|417p . we find that, for V >h, the logarithmic 
terms at ft, = and h = V are essentially proportional to 
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the longitudinal exchange coupling . This means that 
even in the limit — > 0, these logarithmic contributions 
survive, which is an interesting nonequilibrium effect. 



C. Current 

The stationary current can be calculated from 
Eq. (j23T]) . 

{py^ = r^(o) + 2Afr^^(o) , (420) 

where the stationary magnetization M follows from 
Eq. ppg]) and the components r^(0) and £^(0) of the 
current kernel are listed in Eqs. ([335)) . ((5M)) . (|55g)) . 
and (I368D . For the special case (|376p of two reservoirs 



and the relation (|379p between the exchange couplings, 
we obtain 



T'^{Q) = llv[{j:^f + 2{jt^f} 



-pis: 



(0) 



•■'a '^ndj 
2 



(421) 



(422) 



where we have used cj^^, = —^a'ySa^-a'- 

For V <^h, the logarithmic terms '-^ C2{h), ^2{V — h) 
can be neglected and the magnetization is given by 
M = -\. According to Eq. (|420| . this leads to (/^)^* = 
r^(0) — r!^^(0), and we see that the logarithmic terms 
~ Ci{V) cancel. This gives the result of elastic cotun- 
neling 



(423) 



but with renormalized exchange couplings. 

Ai V ~ h, the differential conductance — ^ has 
an interesting feature, see Fig. [T^ For voltages slightly 
below the renormalized magnetic field {V < h — V <^ 
h), we get in units of Go — e'^/h (as before, we list the 
leading order and logarithmic terms in next-to-leading 
order) 



nd) 



C2{V-h) . 

(424) 

For voltages slightly above the renormalized magnetic 
field {V >h,V -h-€.h), we get 



{JI + Jk) 



2J, 



2 7^ 



C2{V~h) 



(425) 



Two interesting features happen at V ~ h: There is a 
jump in the leading order term in the conductance (if 
7^ 0) that is due to inelastic cotunneling which sets in 
at this voltage. The jump is superposed by a logarithmic 
term which becomes largest for V — h. The experimen- 
tally accessible parameters characterizing the line shape 
are given by the position, broadening and height of the 
resonance. The position is approximately at V = h [up 
to unimportant terms of the 0(r2)], and the broadening 
of the left side of the resonance is given by Using 
Eqs. p95p and p94p . we obtain at resonance 



J, 



nd; 



(426) 
(427) 



The right side of the resonance has no characteristic 
broadening: asymptotically it reveals the logarithmic 
voltage dependence of the exchange couplings. The value 
of the conductance at resonance can be calculated in good 
approximation by taking the extrapolation of Eq. (I425P 
at V = h, i.e., by replacing C2{V — h) In^r-r- — in 



- 2ll/ = h 



Eq. (1123 • Together with Eq. this reveals the 

Kondo-induced logarithmic enhancement ^ J'^ In J of the 
conductance at resonance superposed on the enhance- 
ment from inelastic cotunneling. The result for the max- 
imal conductance can be simplified by using the relation 
(|379p . which gives 

(G^/Go) max 

■[{pf + 2{l + 2xLXR){j^f + 



+ 2{3 + AXLXR).r{J^flll-:r- 



V 



, (428) 



with 



ilv^h = {^^LXR {rf + [J^f] . (429) 



This result for the value of the maximum conductance at 
resonance has first been noted in Ref. 19 for the case of 
the isotropic Kondo model and x^ = xr = i. 

Our results for the differential conductance agree with 
those of Refs. [lO and [H. However, note that the pre- 
cise line shape at resonance can only be obtained if one 
includes the cutoff scales T1/2 into the logarithms and 
uses the correct smearing p83p of the sign function. In 
the numerical plots, we have always used these smeared 
forms and we calculated h and T1/2 self-consistenly from 
Eqs. (|393p - (|395p (the same procedure has been used for 
the susceptibility). 



D. Results for the anisotropic Kondo model 

In this section, we discuss how the results for the dif- 
ferential conductance, the susceptibility, the relaxation 



39 



0.4 




I 1 ' ' 1 

0.5 1 1.5 2 

V/ho 

FIG. 12; The differential conductance G = as function of 
bias voltage at ho = 10~^_D for the isotropic Kondo model 
with Tk = W'^D. 



and dephasing rates, and the renormalized g factor are 
affected if we consider the symmetric, but anisotropic 
Kondo model with J^^, = ^ ^ Jaa'- recently 
proposed,— such models can, e.g., be experimentally re- 
alized by studying quantum transport through single 
molecular magnets, where the transverse exchange cou- 
pling is induced by magnetic quantum tunnehng terms 
due to transverse anisotropics. To investigate this, we 
vary and such that the Kondo temperature Tk 
remains constant and only the other invariant of the RG 
equations, i.e., = (J^)^ — (J^)^ is changed, see (|270p . 
If is positive, is larger and is smaller than the 
isotropic coupling in the case = 0. 
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FIG. 13: The differential conductance G = -r^ as function of 

dV 

bias voltage at ho = 10~^-D for the isotropic Kondo model 
(solid line) and the anisotropic Kondo model with two dif- 
ferent values of = (J^)^ — (J^)^ (dashed and dash-dotted 
lines). The Kondo temperature Tk ~ 10^'' D is the same in 
all cases. 



The differential conductance is shown for three differ- 
ent values of in Fig. 1131 For growing c^, i.e., increas- 
ing and decreasing J^, we can observe a number of 
different effects compared to the isotropic case: First, 
the conductance for y ^ ft, which is proportional to 
("^nd)^ according to Eq. (j423p is increased. Second, the 
step at V — h that is due to inelastic cotunneling de- 
creases its height which is proportional to ( Jj^)^ in lead- 
ing order according to Eqs. ()424|) and (|425|) . Third, also 
the logarithmic enhancement at y = ft is reduced for 
increasing c^. The prefactors of the logarithmic terms 
in Eqs. (|424p and (|425p are proportional to J^(J^)^, a 
quantity which decreases for increasing c^. Finally, the 
position of the step and the logarithmic enhancement 
shifts to the right, i.e., h is increased. The reason is that 
the leading contribution to ft — fto is negative and propor- 
tional to — {Ja)o, i-e., the renormalization of the cou- 
pling J^, according to Eq. (|392p . This renormalization 
is in leading order determined by {J^d)^^ see Eq. (|266p . 
An increase in which is connected to a decrease in 
therefore leads to a weaker renormalization of and 
hence to a larger value of h. 

Figure [T3] shows the magnetic susceptibility x = ; 

multiplied with — (^"i^'') ^ for different values of the 
anisotropy. For larger c^, i.e., increased J^, the logarith- 
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FIG. 14: The susceptibility x, muhiplied with -i^^, as 
function of the magnetic field at V — 10~*Z) for the isotropic 
Kondo model (solid line) and the anisotropic Kondo model 
with two different values of = (J^)^ — (J^)^ (dashed and 
dash-dotted lines). The Kondo temperature Tk ~ 10~*-D is 
the same in all cases. 

mic enhancements become larger, see Eq. (j419|) . This 
is an effect already mentioned at the end of Sec. IIVBI 
The susceptibility depends only on the longitudinal cou- 
pling and the logarithmic contributions survive even 
in the limit 0. Especially for molecular magnets 

with a very small quantum tunneling term, this might be 
an interesting possibility to observe logarithmic contribu- 
tions in the susceptibility although the Kondo tempera- 
ture might be quite small. Furthermore, the resonance 
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at the point where h = V shifts sHghtly, this is due to the 
dependence of h on the anisotropy which was discussed 
above. 

The derivatives of the rates Fi and r2 with respect to 
the magnetic field ho are shown in Figs. [15] and [TBI re- 
spectively. Increasing c^, i.e., decreasing J-"-, leads to a 
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FIG. 15: The rate Fi, derived with respect to the magnetic 
field, as function of the magnetic field at V — 10~*D for the 
isotropic Kondo model (solid line) and the anisotropic Kondo 
model with two different values of = {J")^ — (J^)^ (dashed 
and dash-dotted lines) . The Kondo temperature Tk- — lO^^D 
is the same in all cases. 

decrease in the leading-order contributions to the deriva- 
tives of the rates which are proportional to (J,^)^, see 
Eqs. (|396p and (|397p . The dependence of the logarithmic 
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FIG. 16: The rate F2, derived with respect to the magnetic 
field, as function of the magnetic field at V — 10~*D for the 
isotropic Kondo model (solid line) and the anisotropic Kondo 
model with two different values of = {J")^ — (J^)^ (dashed 
and dash-dotted lines). The Kondo temperature Tk ~ lO'^D 
is the same in all cases. 

effects at h = and h = V on the anisotropy is dif- 
ferent for the derivatives of the two different rates: The 



features in dVi /dho become less pronounced whereas the 
resonances in dT2/dhQ become sharper. The reason is as 
follows: all logarithmic terms in next to leading order in 
the rates are proportional to J^{J^^)'^ which decreases 
for larger anisotropy, but the logarithms are broadened 
by r2 in the case of dTi/dho [see Eq. (I396p ] and by Fi 
in the case of dT2/dhQ [see Eq. (|397p ]. The fact that Fi, 
which is proportional to {J^Y' in leading order accord- 
ing to Eq. (|393p . decreases for larger c? , and F2, which 
contains a term oc ( J^)^, see Eq. (|394p . increases with c^, 
explains the observed behavior. 

Finally, we discuss how the renormalized g factor is af- 
fected by the anisotropy, see Fig.[T71 For growing c^, the 
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FIG. 17: The renormalized g factor g = 2dh/dho as function 
of the magnetic field at V = 10~*D for the isotropic Kondo 
model (solid line) and the anisotropic Kondo model with two 
different values of = (J^)^ — (J^)^ (dashed and dash-dotted 
lines). The Kondo temperature Tk ~ 10~^D is the same in 
all cases. 

deviation of g from the unrenormalized value 2 becomes 
smaller. The reason has already been discussed above: 
Like the renormalized field h, the g factor (|398p is in first 
order given by the renormalization of the coupling 
which is determined by ( J^)^ and is therefore decreased 
for positive c^. Also the prefactors of the logarithmic 
terms depend on (J^)^ and therefore become smaller for 
growing c^. 



V. SUMMARY 

We have presented a systematic and analytic weak- 
coupling analysis in two-loop for a generic quantum dot 
coupled to reservoirs via spin or orbital fluctuations. We 
have used a nonequilibrium RG formalism in Liouville 
space (or, equivalently, on the Keldysh contour), which 
is a natural and formally exact generalization of conven- 
tional poor man scaling methods. The essential difference 
between the RTRG-FS method and other RG formalisms 
in nonequilibriuniiSiS£'2l, is the fact that RTRG-FS ana- 
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lyzes the dynamics of the reduced dot density matrix in 
Laplace space by calculating an effective dot Liouvillian. 
In each step of the RG approach an effective kinetic equa- 
tion is obtained. The RG equations are coupled differ- 
ential equations for the kernel determining this kinetic 
equation and for the vertices in Liouville space. The 
structure in Liouville space is essential for the genera- 
tion of dissipation and to obtain the transport relaxation 
and dephasing rates together with the transport oscil- 
lation frequencies of the various exponentially damped 
modes (rather than quantities of quantum decay the- 
ory, which describe the decay of a discrete dot state into 
the reservoir continuum). Technically, it is known from 
the Keldysh formalism^ that dissipation (or the collision 
term of a quantum Boltzmann equation) is generated by 
diagrams connecting the upper with the lower part of 
the Keldysh contour. Therefore, to get the correct trans- 
port rates and the precise position of resonances from 
the oscillation frequencies, it is essential to include these 
contributions in the theory and, at the same time, com- 
bine it with RG methods to get rid of logarithmic diver- 
gencies. We note that the RTRG-FS formalism is not 
based on the usual Keldysh formalism as, e.g., methods 
based on slave particles^ or nonequilibrium functional 
RG methods.— Since the unperturbed part of the Hamil- 
tonian contains strong interactions, the RTRG-FS inte- 
grates out only the noninteracting reservoirs but takes 
the correlation effects on the quantum dot exactly into 
account. This leads to a natural formalism in terms of 
Liouville operators and vertices in Liouville space, which 
we have combined with quantum field theoretical meth- 
ods of diagrammatic representations and RG. In this 
sense, we have combined the advantages of Liouville oper- 
ator techniques (providing very compact notations for the 
time evolution of the reduced density matrix in Laplace 
space) , diagrammatic methods (to identify and calculate 
irreducible blocks leading to the oscillation frequencies 
and the relaxation/dephasing rates), and RG methods 
(reorganizing the perturbation theory in such a way that 
logarithmic divergencies are absent). In particular, it is 
possible to prove generically in all orders of perturbation 
theory and within all truncation schemes that the cut- 
off scales for the vertices are the physical relaxation and 
dephasing rates of the exponentially decaying modes of 
the reduced density matrix of the dot. This provides a 
reliable calculation of the line shape at resonances. In 
addition, one obtains the kernel of the kinetic equation 
in Laplace space from which the time evolution into the 
stationary state can be calculated. It is important to 
note that the possibility to calculate the time evolution 
is closely related to the fact that we obtain a systematic 
theory for the line shape at resonances, since the oscilla- 
tion frequencies and the relaxation/dephasing rates gov- 
erning the exponential decay of the time evolution are 
the same scales which determine the resonance positions 
and the cutoff of logarithmic divergencies at resonances 
for stationary quantities such as, e.g., the magnetization 
or the conductance. 



In cornparism to previous formulations of RTRG 
(Refs. [25l - [27l) (where the physics of cutoff scales by relax- 
ation/dephasing modes was also included), the recently 
proposed RTRG-FS method^- is a formulation in pure 
Matsubara frequency space and avoids the need of the 
Keldysh indices after a first discrete RG step where the 
symmetric part of the Fermi distribution is integrated 
out. The main advantage of this development is the 
possibility to prove that the perturbation theory in the 
renormalized vertices on the r.h.s. of the RG equation is a 
well defined series in the weak-coupling regime. Essential 
for this proof is the fact that the imaginary part of the 
denominators of all resolvents consists of a sum of strictly 
positive Matsubara frequencies (bounded from above by 
the cutoff), the cutoff scale A, and some positive relax- 
ation or dephasing rate F^. As a consequence, it is pos- 
sible to solve the RG equations analytically in the weak- 
coupling regime by a systematic perturbative expansion 
in the renormalized couplings. In this paper, we have 
described this procedure for a generic quantum dot in 
the Coulomb blockade regime where spin and/or orbital 
fiuctuations dominate the transport properties. We have 
included all two-loop contributions to calculate physical 
quantities up to the first logarithmic contribution at res- 
onances. We have shown that the latter are well defined 
after renormalization and can be treated within pertur- 
bation theory in the renormalized couplings. 

We have applied the formalism to the anisotropic 
Kondo model at finite magnetic field and finite bias volt- 
age. We calculated stationary properties (magnetic sus- 
ceptibility and conductance) and the parameters charac- 
terizing the exponentially damped modes of the time evo- 
lution (renormalized magnetic field, spin relaxation and 
dephasing rate) . In comparism to the case of zero mag- 
netic field (which was treated within flow equation meth- 
ods already in Ref. Q for the isotropic Kondo model), 
the central issue is the calculation of logarithmically en- 
hanced contributions at the resonance positions h ^ V 
and h <^V. The one- loop leading order theory of Ref.[20l 
was able to obtain these logarithmic contributions but it 
was not shown how to include microscopically the cut- 
off scales from relaxation and dephasing modes into the 
RG formalism. Furthermore, the renormalization of the 
magnetic field was not taken into account and the cal- 
culation was restricted to the stationary values of the 
magnetic susceptibility and the conductance (the spin 
relaxation and dephasing rates were calculated only per- 
turbatively in the bare coupling and in the absence of 
a magnetic field in Ref. ^2^. In this paper, we have in- 
cluded the renormalization of the magnetic field in 0{J) 
and 0{J^ In) together with the spin relaxation and de- 
phasing rate in 0{J^) and 0(J'^ In). This provides a 
consistent theory for the resonance position and the line 
shape at resonance, and answers the question which rate 
occurs in the various logarithmic terms. For stationary 
quantities, the obvious physical result that the spin de- 
phasing rate cuts off the logarithmic contributions is ob- 
tained. For the rates and the renormalized field itself. 



42 



it turns out that the spin relaxation rate determines the 
resonant Hne shape for the spin dephasing rate and the 
g factor, whereas the spin dephasing rate controls the 
resonant line shape for the spin relaxation rate. In ad- 
dition, it turns out that the logarithmic enhancements 
or suppressions at the resonance h — V are quite dif- 
ferent for the various physical quantities regarding the 
symmetry ioi h > V and h < V. The g factor shows a 
symmetric suppression, the conductance an asymmetric 
enhancement, the suspectibility shows only an enhance- 
ment for h < V, the derivative of the spin relax- 
ation rate only an enhancement for V > h, and the cor- 
responding quantity for the spin dephasing rate only a 
suppression for V < h. We proposed an experimental 
setup with a weakly coupled probe lead to measure the 
logarithmic suppression of the renormalized g factor as 
function of voltage V or magnetic field ho. This has to 
be contrasted to the usual measurement of the g factor 
at resonance V = h from the position of the cotunneling- 
and Kondo-enhanced conductance peak, which provides 
the g factor only for one special value of the voltage. In 
the anisotropic case, we find the obvious result that all 
logarithmic contributions involving the transverse cou- 
pling are suppressed for increasing longitudinal coupling 
since the renormalization of the transverse coupling is 
reduced. However, for the magnetic susceptibility this 
does not hold since the transverse coupling cancels out 
in the prefactor of the logarithmic terms. Therefore, the 
logarithmic enhancements become larger for the suscep- 
tibility if one increases the longitudinal coupling. This 
result might be helpful to find Kondo physics in trans- 
port through single molecular magnets, where the trans- 
verse coupling can be smaller than the longitudinal one. 
Concerning the line shape at resonance when increasing 
the longitudinal coupling, it turns out that the logarith- 
mic features become less pronounced if the spin dephas- 
ing rate determines the line shape, whereas it becomes 
sharper when the spin relaxation rate is involved. The 
reason is that the spin dephasing rate involves a voltage- 
induced term ~ V(J^)'^, which increases with increasing 
longitudinal coupling, whereas the spin relaxation rate 
involves only the transverse coupling in order O(J^). 



VI. LIST OF SYMBOLS 

In this section we provide a complete list of all symbols 
used in this paper. We present two tables, one for the 
symbols used in Sec. HI] of the generic case, and one for 
Sec, mil where the method is applied to the Kondo model. 

The table for the symbols used in Sec. [Til is given by 



Symbol 


What the symbol means 


Rcf. 


H 


Hamiltonian 


m 




reservoir Hamiltonian 


111 


Hs 


dot Hamiltonian 


III 


V 


coupling dot ^ reservoir 


m 


Ho 


Hrea + Hs 


m 




frequency 





V 

a 
a 
1 
1 

<5ii' 
ai 

9ii' 

pres 

r„ 

D 

P 
'ii' 

Ps{t) 

[A,B]± 

L 

Lj-y 

Ps{E) 

iniE) 

Trs 

(0) 



Lv 



^11' 

J? 

S 

Np 

X, 

7ii' 

ni...„ 

-El.. .11 

iOl...n 
fj-i 
LUi 

LI 
E(£) 



re// 



(nil. 
pf 

n{z) 

A,(2) 
Zi 

ai 

r 

Gil' 
Gil' 

s I a 

7/ 

r a 

Ga 
11' 

A 
Ao 

At„ 



creation/annihilation index 

index for a,a, . . . 

spin index 

reservoir index 

index for r/vu} 

index for —r]uu) 

field operator 
coupling vertex 
reservoir density matrix 
temperature of reservoir a 
chemical potential of reservoir a 
Fermi function 
reservoir band width 

current operator for reservoir 7 
current vertex 

-|(r?<5a7 -I- rj'Sa'j) 
reduced dot density matrix 
AB±BA 
[H,-]- 

Laplace transform of ps(t) 
Laplace transform of {I''){t) 
trace over the reservoir 
trace over the dot 

Hs,-]~ 
Vr]- 

Keldysh index 

vertex in Liouville space 

field operator in Liouville space 

symmetry factor for equivalent diagrams 

number of permutations of field operators 

energy varia bles crossed by vertical line 

contraction Jf J^, 



L 



{El...n + Ul. 
{El...„ + LUl. 

En — 
,=i^» 

riiLOi 

efi'ective Liouville operator 



Ls{E-i...n + UO\...n)) 



Llff{E)- 



r{0) 



re/// 



current kernel 

current vertex in Liouville space 
stationary dot density matrix 
(z-P/fiz))''^ 
eigenvalues of (z) 
right eigenvector of L^^^ (z) 
eigenvec. with zero eigenval. of P/ 
left eigenvector of L"/^ (z) 
poles of n(2): Zi = \i{zi) 
residua of n{z): (1 - 
Il{z)^ Z/{z~Ls) 

Zi — hi iTi 

short hand notation for Fi = T 

symmetric/antisymmetric part of 7^^^, 
efi'ective Liouvillian including 7^ 
efi'ective kernel including 71 
effective vertex including 7f 
RG cutoff parameter 
initial cutoff 

Matsubara freq. cj" lying closest to 



(131 
(131 

(El 
(131 
(Hi 
(1^1 
(1^1 
(131 
(131 

(Bi 
(151 
(Ui 
(1^1 
(Hi 
(171 
(gQi 
(|li 
(221 
(ll 
(251 
(2|l 

(2^1 
(|gl 
(271 
(261 
(271 

(2' 

(281 
(291 

(m 

(291 

(211 
(311 
(311 

(Ml 

(Ml 
(811 
(^1 

(agi 

(371 
(2Zl 

(^1 

(401 

(Hi 
(Ml 
(|Zi 
(Ml 
(Ml 
(SZi 
(^1 
(|li 
(^1 
(Ml 
(561 

dMl) 

(Hi 
(731 
(Z|i 
(751 

(gi 
IB3I1 
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Ac 
I 

Ls{E,uj) 

Gl2(...) 

Ja 

Jc 

Jo 

Ka{,z) 
KAiz) 
Fa{z) 
Fa{z) 
h{z) 

■-^12 



ri)c 



M2a/b) 
^{2<ii/2) 

^(26/c) 
r(3a/i)) 

Ai...„ 
A 



max{|i5|, \hi\} 
In^ 

cutofT-dependent Fermi function 
Matsubara frequencies for reservoir a 
smeared ^-function 
cutoff-dependent antisymmetric part 
Ls{E + iuj) 

{E + iLo-Ls{E,Lu))-^ 
order of Gn/ ~ 0{Ja) 

Ja=Ao 

^a(^) - f 

7?aW = ^^aI^) 
FaW- f ln2 



0(J") of Ls(£,<^) 

0(J") of Gl2(S,<^,C^l,C^2) 



lI^' (£, tj) = L'g' -{E + iuj)Z 
Ip{E,co)\A=A^ 

G]^\E,ljJ,LUi,LLl2)\A=Aa 

^(2) _ M2a) M2b) 
'-'11' - '-'11' + '-'11' 
M2a) _ ■M2ai) M2a2) 
'-'11' - *'-'ll' +'-'11' 

decomposition of L''? {E, lu) 

decomposition of L? {E, u)) 

decomposition of Lg {E, u)) 

Ei,,,n + iuj — Ls 
short hand notation for Ai, 



(1) 



7(1) 



A 



( TTTj) 
266)) 

(Si 
(761 
T7 



(Ml 
(911 



(fT23)) 
((T82|) 
(fT25|) 

(nssi 

(fT33)) 
(fT34)) 
(fT29l) 



(U36|) 

mil) 

(fT56)) 

nMi) 

(p76l) 

unsi) 



The table for the symbols used in Sec. lIIII is given by 



Symbol 



ho 
h 
h 
fi 

f2 

s 

P , 

^ act' 

j./± 

jz/± 
jz/± 



Tk 
T'k 

V 

L± 

L" 
L" 
L"^ 
Li 
L2 
Ls 
i± 
Li 
Ll 



What the symbol means 



bare magnetic field 
(1 - - iJ^)o))ho 
renormalized magnetic field 
spin relaxation rate 
spin dephasing rate 
spin operator 

bare exchange couplings (i — x, y, z, ±) 
renormalized exchange couplings 



aa' 



--/z/± 

-y' O 



O 

act 

•'LR — •'RL 

factor for coupling to reservoir a 
Kondo temperature 
voltage 

L+A = SA,L-A 

I- 1 + L+-L- 

II- L+-L- 
I • 1 + 2L+'L-' 

Wr+ L- -2iL 



-AS 



2 



- \ L - 

Ll± 
LlT 



+ L-) 
- L- + 2iL^ 
= 1,2. 
+ L+ 



{LtL- 
{LlL- 



X L-) 
3,± 
LI) 



+ L+^Ll) 



Ref. 



(1991 



( 392 1 



( 335 1 



( 334 1 



( 200 1 



( 266 1 



( 268 1 



( 377 1 



( 265 1 



( 265 1 



335 



200 



377 
203 



( 270 1 



(3161 



( 270 1 



( 204 1 



( 207 1 



( 208 1 



( 208 1 



( 208 1 



( 208 1 



( 209 1 



( 209 1 



( 209 1 



(2101 



h{E,co) 

r(E,Lj) 

r\{E,co) 

ME,oj) 

Pr{E,w) 

\^r{E,i^)) 

{HE,io)\ 
zo 

Zl 

z± 

<7± 

LI 
LI 

G''=(...) 
G^=(...)^ 

{L-r" 

C^[X) 
rp 

^aa' 

sign,(a::) 

\x\ ■ 

e^{x) 

PUl 

M 

X 

g 

X 
Go 



term h{E,Lo)L^ in Ls{E,lo) 

term —iV''{E,uj)U in Ls{E,uj); i = a,c,3z 

term iV\{E,ixi)V in T.^{E,uj); i^b,lz 
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APPENDIX A: SECOND-ORDER RG 

In this appendix, we prove that the effect of the RG 
terms (|149p and (|150p is the replacement F\{z) — * F^{z) 
in Eq. (fT5^ . with Fj^{z) defined by Eq. (fT%)) . 

First we consider the term (|150p . The terms of 0(-^ J'^) 
can be extracted by using the expansion (jl5ip for the 
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integral over the resolvent and considering only the terms i.e., the terms cancel each other, 
of 0(1) and 0(4) in this expansion, 



I I dw2n(£;i2,A + a> + W2) »ln(2) + — zi2, (Al) 



where we have defined 



In contrast, Eqs. (jA5p and (|A7|) are not independent of 
the integration variable 072- When inserted in Eq. ()149|) . 
we need the integral (zy is replaced by its eigenvalue in 
this equation) 



zi2 = {E 12 + iuj - L^g'') 



(A2) 



Inserting this in Eq. pSOp for the two independent in- 
tegrals over the resolvents and collecting the terms in 
0(^ J3) leads to the result 

(fT50D 

--^^ |g'i2''zi2G'23 Gg^ + Gi^G^g zisGgj^l . (A3) 

Next we consider the term (|149p and use for G'^^' the 
result (|142p . which gives 

^(2b) 



G^f^(£;i2,A + co + u;2,-^2,-A) 



rW,„ 2A-ZZ13 -(1) 

•-^23 X 31 



and 



(A4) 



.G«lnA±^^^|^G^|), (A5) 



duJ2 -; In 

A + UJ2 - izi2 



K + UJ2- iZ23 

A 



1 , 2.„^ ln(2) i , 

^ 2^^" (2)-*^^i2+2a(^i2-Z23), 
where we have expanded up to linear order in This 
gives 



(IT491) from ([All and (jA7| 



ln(2) 
2A 



\<^12 ^12<^13 <^32 + <^23 <^31 ^12*^21 / 

1 //^(l)^ , ^(1)^(1)^ (=,{\) 

-— |Gri2 ^12tTj3 LT32 + Lr23 '-'31 ^12^21 



"'-^12 '-'13 ^23Lf; 



32 '-'23 



,;,)z23G«G«} (AlO) 



Using the antisymmetry property ([97)) . we see that the 
second term is zero and the first term agrees with 
Eq. (jA3p . leading to the final result 



Glf (£;,a.,A,c.2) = G«ln^A:p^G(^,) (A6) 

Equations (jA4p and (jA6P are independent of the in- 
tegration variable 072, and we can replace the logarithm 

by 

ln^^^^«ln(2)-^zi3. (A8) 

Inserting Eqs. (|A4| and (|A6| in (fTiO]) . using Eq. ((AT|) . 
and collecting all terms of 0(-j J'^), we find 

(fT49| from ([Ail and (|A6| ^ 

M^r^(i)^ ^(1)^(1) ^(1)^(1)^ ^(1) 

1 '-'12 ^12'-^23 31 ~ 12 23 ^IS'-'gj 



2A 



G 



Zi3 0^2'' ^^l'' + ^13 ^32^^12^21''} - 0'(A9) 



ln(2) 
A 



|Gi2 (-£^12 -^lU) - Lg jG23 Ggj 



■ 6^2^43 (^^13 + *^- 4°^ )G^i^} • (All) 



Finally, using the leading-order RG equations (|115p or 
(I116p . we see that Eq. (jAlip can be written in the form 



(HHl) + (HSOl) - 

ln(2) d 



2 dA 



{Gl^2^(£;i2+i^-4"^)G^l^} , (A12) 



i.e., a similiar term to Eq. (|152p results, which proves the 
replacement (I155p . 
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